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ABSTRACT 


This  thesis  provides  a  crosstalk  analysis  of  optical  chip  interconnects  via  single- 
mode  waveguides  with  synchronous  transmission  and  asynchronous  transmission.  This 
crosstalk  noodel  is  general  and  can  be  used  for  any  type  of  waveguide  network.  Three 
cases  of  laser  sources  will  be  considered:  (1)  each  channel  operates  with  an  independent 
laser  sources,  (2)  all  laser  sources  have  the  same  mean  wavelength  but  with  different 
phase  noise  processes,  and  (3)  all  laser  sources  are  identical  with  the  exception  of  the 
initial  phases.  The  analysis  takes  into  account  the  coupling-induced  crosstalks  between 
adjacent  waveguides,  the  laser  linewidth.  the  shot  noise,  the  dark  current  generated  by 
the  photodiode,  and  the  post-detection  thermal  noise.  Bit  error  probabilities  versus 
received  peak  powers  are  presented  together  with  power  penalties. 
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I.  INTRODUCTION 


This  thesis  considers  the  problem  of  crosstalk  impact  in  an  optical  interconnect 
system  using  a  single-mode  waveguide  network  [Ref.  1-Ref.  8].  An  optical  interconnect 
consists  of  two  or  more  terminal  nodes  which  are  interconnected  hv  a  single  optical  path 
or  a  high-density  parallel  network  which  preserves  the  parallel  nature  of  the  data 
generated  at  the  nodes  [Ref.  9].  Such  interconnect  systems  would  consist  of  chip-to-chip 
or  boaid-to-board  interconnections  [Ref.  8],  The  use  of  waveguides  provides  a  potential 
of  integrating  the  entire  system  of  transmitter  and  receiver  on  the  same  substrate  [Ref. 
10]. 


XMTR 


RCVR 


Figure  1:  Schematic  diagram  of  an  optical  interconnection  transmission  system 

architecture  [  Ref.  9]. 


A.  BACKGROUND 

The  success  of  optical  communication  has  accelerated  research  on  high  capacity  data 
handling  systems.  It  is  expected  that  the  monolithic  integration  of  optical  and  electronic 
components  on  the  same  chip  will  ultimately  lead  to  ultrahigh-speed,  high  sensitivity, 
compactness,  reliability,  low  cost,  as  well  as  passive  and  active  integrated  optic 
components  [Ref.  11]. 

Integrated  optics  refers  to  the  implementation  of  various  functions  with  light  such  as 
modulation,  switching,  generation,  and  detection  in  an  optical  guided  wave  structure 
formed  on  a  substrate.  Use  of  the  term  "integrated"  implies  implementation  of  many  of 
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these  functions  on  the  same  substrate.  Integrated  optical  devices  are  characterized  by 
many  advantages  associated  with  lightwave  technology  in  general,  namely,  larger 
information  capacities  than  with  electronic  technology  are  possible,  electromagnetic 
interference  is  not  present,  and  parallel  processing  is  possible.  Many  of  the  advantages  of 
electronic  integrated  circuits,  including  the  potential  for  fabrication  economy  and 
reliability  of  devices  combined  onto  one  substrate  apply  to  integrated  optics.  Vibration 
problems  associated  with  bulk  optical  experiments  are  eliminated  when  devices  can  be 
integrated  onto  one  substrate  space.  [Ref.  12] 

Qosely  associated  with  integrated  optics  is  the  field  of  integrated  optoelectronics 
which  encompasses  device  structures  referred  to  as  optoelectronic  integrated  circuits 
(OEICs)  [Ref.  12].  OEICs  represent  a  device  technology  with  potential  to  meet  a  broad 
range  of  future  telecommunication  and  computing  system  needs.  Optoelectronic 
integrated  circuits  are  circuits  that  monolithically  integrate  optical  and  electrical 
components  on  a  single  semiconductor  chip  [Ref.  10].  Compound  semiconductor 
materials  are  used  to  form  optoelectronic  integrated  circuits  because  photonic  devices 
such  as  semiconductor  lasers,  detectors,  high-speed  electro-optic  modulators  and 
switches,  as  well  as  quantum  well  waveguide  devices,  and  high-speed  electronic  devices 
can  be  formed  with  the  same  material  alloys  [Ref.  12]. 

High  performance  is  now  being  achieved  in  devices  using  both  GaAs  and  InP 
material  systems  [Ref.  9].  It  is  in  the  combination  of  photonic  components  and  electronic 
circuitry  that  the  OEIC  gains  a  usefulness  over  and  above  that  which  can  be  obtained  by 
placing  non-monolithic  circuits  together  in  a  package.  For  example,  it  has  been  assened 
that  only  through  monolithic  integration  can  one  fabricate  extremely  high-bandwidth 
transmitters  or  high  sensitivity  receivers  [Ref.  9]. 

OEICs  currently  fall  into  three  categories:  interconnects,  communications,  and 
computing  and  signal  processing.  As  mentioned  before,  an  optical  interconnect  (Fig.  1) 
consists  of  two  or  more  terminal  nodes  which  are  interconnected  by  a  single,  high-band 
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optical  path,  or  alternately,  a  high-density  parallel  optical  network  which  might  preserve 
the  parallel  nature  of  the  data  generated  at  the  terminals.  Such  an  architecture  is  useful 
for  interconnecting  two  computer  mainfnunes  or,  on  a  smaller  scale,  would  consist  of 
board-to-board  or  even  chip-to-chip  interconnections.  The  desirability  of  such  a  system 
is  its  potential  for  the  rapid  transfer  of  parallel  data  from  one  system  to  another  without 
the  concomitant  problems  of  electromagnetic  interference  ap  j  signal  dispersion 
commonly  observed  in  a  high-bit-rate  electrical  interconnects.  [Ref.  9] 

Fci  example,  as  computing  environments  evolve  toward  the  type  of  distributed 
network  in  Fig.  2  with  data  processing  and  database  sharing  among  remote  locations, 
overall  operation  efficiency  increasingly  relies  on  the  efficiency  of  interconnection  links. 
Optical  interconnects  are  a  new  approach  that  could  be  used  to  achieve  high-bandwidth 
low-loss  interconnects  for  these  applications.  [Ref.  9] 

Figure  3  illustrates  the  functional  charac  .eristics  of  a  GaAs  computer  interface  chip  . 
For  insertion  into  a  practical  computer  system,  this  circuit,  for  any  OEIC  component,  has 
been  designed  to  meet  the  following  requirements  [Ref.  10]: 

1)  The  circuits  must  be  high  speed,  capable  of  multi-gigabits  per  second,  with 
relatively  high  complexity  (lOK  transistors/chip). 

2)  They  must  be  capable  of  operation  in  noisy  environments;  Pj,  <  lO’^^  in  the 
presence  of  noise  levels  up  to  100  mV. 

3)  They  must  have  high  reliability,  i.e.,  a  failure  rate  <0.01%  per  kilohours  for  the 
entire  link,  for  over  10^  h  at  50®  C. 

4)  They  must  have  redundancy  in  critical  paths  with  transmission  error  correction 
capacity. 

5)  They  must  be  compatible  with  existing  computer  interface  technologies, 
packaging,  power  supply  requirements,  and  silicon  interface. 
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Figure  2:  Distributed  computer  network  [Ref.  10]. 


The  advantages  of  OEICs  are  currently  being  exploited  by  the  Department  of  the 
Defense.  OEICs  are  being  developed  for  such  applications  as  neural  networks,  high 
speed  signal  processing,  high  speed  communications,  and  advanced  antenna  systems.  For 
example,  the  development  of  efficient  microwave  signal  distribution  by  optical  fiber 
creates  the  possibility  of  achieving  unconventional  antennas  such  as  "sman  skins".  In 
this  example,  individual  radiated  elements  would  conform  to  the  contour  of  an  aircraft 
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and  the  "antenna's"  directional  propenies  would  be  determined  by  control  of  phase 
emission  from  different  parts  of  the  aircraft.  Consequently,  the  DOD  funds  on-going 
programs  developing  high  frequency  components,  such  as  high  speed  GaAs  circuitry. 
This  program  is  funded  for  seven  years  at  a  level  of  $500  million.  [Ref.  10] 


Figure  3:  OEIC  Functions  [Ref.  10] 


B.  APPROACH 

This  thesis  considers  a  mathematical  approach  that  includes  moTe  than  two 
intnfoing  channels.  The  modulation  format  used  in  this  analysis  is  on-off  keying  (OOK) 
wdth  direct  detection.  Two,  three,  and  four  channels  were  modeled  using  synchronous  bit 
transmission.  For  comparison,  two  adjacent  channels  were  modeled  using  asynchronous 
bit  transmission.  The  emphasis  of  this  thesis  will  be  on  the  impact  of  coupling  induced 
crosstalk  between  adjacent  waveguides  in  a  waveguide  network.  Furthermore,  additional 
parameters  that  will  affect  the  performance  of  the  interconnect  system  are  the  laser 


linewidth,  the  shot  noise  and  the  daric  current  generated  by  the  photodiode,  and  the  post- 
detection  thermal  noise.  To  be  general,  the  receiver  will  be  considered  to  be  an  integratc- 
and-dump  filter  with  integration  time  T  where  T  is  the  bit  time.  Also,  the  current  spectral 
density  of  the  post-detection  thermal  noise  is  denoted  as  Nq.  The  spectral  density  Nq  can 
be  easily  computed  for  a  given  low  noise  amplifier  type  given  the  effective  noise 
temperature  and  the  matched  resistance  load. 

The  optical  interconnect  system  is  given  in  Fig.  4(a)  for  CX)K  direct  detection.  The 
waveguide  network  consists  of  many  single-mode  waveguides.  It  can  be  a  planar  array 
with  uniform  waveguide  separation  or  any  other  structure  (Fig.  4(b)).  The  envelope  of 
the  output  lightwave  of  a  given  waveguide  is  detected  by  the  photodiode,  which  also 
generates  shot  noise  and  dark  current.  The  output  of  the  photodiode  is  further  corrupted 
by  the  amplifier  thermal  noise.  The  total  signal  plus  crosstalk,  shot  noise,  dark  current, 
and  thermal  noise  is  integrated  over  one  bit  time  T  and  the  resulting  bit  energy  at  the  end 
of  each  integration  time  is  determined  by  the  sheer  to  be  either  bit  one  or  bit  zero.  It  is 
assumed  that  the  waveguide  bandwidth  is  much  larger  than  the  bit  rate  and  the  signal 
spectrum. 

The  resulting  model  for  synchronous  and  asynchronous  transmission  was  validated 
using  MATLAB.  The  resulting  data  from  the  MATLAB  programs  was  graphed  as  a 
series  of  performance  and  optimal  threshold  curves  versus  peak  powers.  These  resulting 
curves  were  then  compared  with  Ref.  8  and  Ref.  13  to  determine  the  credibility  of  the 
results.  Appendices  C  through  F  contain  sample  programs  for  synchronous  and 
asynchronous  transmission. 
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Figure  4(b):  Planar  Waveguides 


C.  THESIS  ORGANIZATION 

This  thesis  is  organized  into  four  chapters.  Chapter  I  provides  a  background  to 
optical  chip  interconnects.  Chapter  II  provides  the  theoretical  framework  for  the 
performance  analysis  and  discusses  the  numerical  results  when  the  bit  transmission  is 
assumed  to  be  synchronous  where  bit  streams  in  all  channels  are  time-aligned.  Chapter 
m  provides  the  theoretical  framework  for  the  performance  analysis  and  discusses  the 
numerical  results  when  the  bit  transmission  is  assumed  to  be  asynchronous  where  bit 
streams  in  all  channels  are  not  time  aligned.  As  in  Ref.  8,  three  cases  of  laser  sources 
were  considered  and  modeled.  In  case  one,  the  system  was  modeled  utilizing  non- 
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coinciding  and  uncorrelated  light  sources  (each  channel  operates  with  an  independent 
laser  source).  In  case  two,  the  system  was  modeled  utilizing  coinciding  and  uncoirelated 
channel  light  sources  (all  laser  sources  have  the  same  mean  wavelength  but  with  different 
phase  noise  processes).  Finally,  in  case  three,  the  system  was  modeled  utilizing 
coinciding  and  correlated  channel  light  sources  (all  laser  sources  are  identical  with  the 
exception  of  the  initial  phases).  Chapter  IV  is  a  summary  of  the  study.  Appendices  A 
and  B  are  derivations  of  the  conditional  mean  square  value  of  Xi  for  synchronous  and 
asynchronous  transmission.  Appendices  C  through  F  contain  sample  MATLAB 
programs.  Appendices  C  and  D  are  sample  programs  for  two  adjacent  channels  with 
synchronous  and  asynchronous  transmission.  Appendices  E  and  F  are  sample  MATLAB 
programs  for  three  and  four  adjacent  channels  with  synchronous  transmission. 
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n.  SYNCHRONOUS  OPTICAL  CHIP  INTERCONNECTS 


A.  ANALYSIS 


For  mathematical  convenience,  the  complex  envelope  notation  of  a  real  signal  is 
adopted.  Thus,  for  a  given  transmitted  bit  bjo  of  a  given  channel  0  whose  laser  phase 
noise  process  is  6o(t) ,  the  signal  at  the  input  of  the  photodiode  is  designated  as 


Si(t)  =  ^bioC^o(‘)  +  0  S  t  S  T 


(1) 


where  A  is  the  OOK  signal  amplitude,  6j;(t)  is  the  laser  phase  noise  process,  coj^  and 
are  the  frequency  spacing  and  the  initial  phase  difference  between  channels  k  and  0, 
respectively.  Also,  bj]^  represents  the  transmitted  bit  of  channel  k,  and  represents  the 
coupling  from  channel  k  to  channel  0.  The  summation  term  in  equation  (1)  thus 
represents  the  crosstalk  from  M  adjacent  channels  into  channel  0.  In  equation  (1),  all 
waveguides  are  assumed  to  have  the  same  attenuation.  For  the  case  of  M  2,  that  is, 
when  only  the  two  nearest  adjacent  channels  are  considered,  equation  (1)  reduces  to  the 
result  in  equation  (2.4)  in  Ref.  8.  In  this  case  =  b/Vi-2B^,  k=l,2  where  B  is  the 
total  power  coupled  from  a  waveguide  into  its  adjacent  waveguide.  In  practice  bj^,  k  = 
0,l,—>i  are  not  necessarily  equal  to  zero  or  one.  Let  r^,  k  ==  0,1,'",M  be  the  extinction 

ratio  of  the  laser  of  channel  k  defined  as  the  ratio  of  the  transmitted  power  of  the  logical 
zero  to  that  of  a  logical  one.  Then  b,^ ’= /(!'*' logical  zero  and 
bik  =  .^1/(1 +  rk)  for  logical  one.  Hereafter,  the  t^'s  are  assumed  to  be  identical  for  all 
channels. 

Let  R  =  neq/hf  be  the  photodiode  responsidvity  [Ref.  14  -  Ref.  15]  where  n^  S  1  is 
the  quantum  efficiency,  q  is  the  electron  charge  (1.6x10-19  h  is  Planck's  constant 
(6.626  X  10-34  j.s),  and  f  the  frequency.  The  output  current  of  the  photodiode  is 
*^l3i(t)l  +  w,(t)+w^(t)  wh^e  w,(t)  is  the  shot  noise  generated  by  the  photodiode  and 
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w^kCt)  is  the  dark  current  noise.  The  output  of  the  photodiode  plus  the  post-detection 
thermal  noise  n(t)  is  integrated  by  the  integrate-and  dump  filter  with  the  normalization 
constant  R  resulting  in  the  following  decision  variable  Yj 


=  iJo  ^ ^ Jo  *dk(>)<lt +ij„’^n(t)dt 

=  Xi  +  W,  +  Wdk  +  N 

where 

Xi  =  Jo|si(t)|'dt 

W,=-^JoW,(t)dt 

Wdk=|;jJwdk(t)dt 

N  =  — fjn(t)dt. 
rJo  ^ 


(2a) 


(2b) 

(2c) 

(2d) 

(2e) 


Since  n(t)  is  a  zero  mean  Gaussian  process  with  spectral  density  No,  the  Gaussian  random 
variable  N  also  has  zero  mean  and  its  variance  oj^  is  given  by 


(3) 


On  the  other  hand,  the  shot  noise  w,vO  is  a  non-stationary  process  since  the  envelope  of 
the  signal  at  the  input  of  the  photodiode,  namely  |si<t)|,  is  time-dependent.  Therefore,  the 

shot  noise  w,(t)  can  be  modeled  as  a  zero  mean  wide-sense  stationary  Gaussian  process 
whose  spectral  density  Wo(bi),  given  a  bit  pattern  bi  =  (bjo,  bii,”',bik),  is  proportional  to 

the  conditional  mean  of  the  squared  envelope  of  the  input  signal.  In  other  words, 
E{lsi(t)l^lbi|  approximates  lsi(t)p  over  a  bit  time  T.  Based  on  this  approximacon,  the 
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shot  noise  spectral  density  can  be  obtained  as  follows  [Ref.  17] 


W„(bi)-<,RE{ls,(l)|^|bi} 

=  iqRA»fb5,  +  £(^bL 

^  V  k»i  y 


(4) 


From  equation  (2c)  and  equation  (4),  the  conditional  variance  of  W,.  given  a  bit  pattern 
bi.is 


(5) 


The  daiic  current  noise  spectral  density  function  is  qhjk  where  Idk  is  the  dark  current.  The 
variance  of  the  dark  current  noise  is  =:TqI^/R^.  The  random  variable  Xi  in 

equation  (2b)  consists  of  the  signal  term,  the  signal-crosstalk  term,  and  the  crosstalk- 
crosstalk  term.  Substituting  equation  (1)  into  equation  (2b),  the  random  variable  becomes 

X,  =  i A^Tbfo  +  i  2CkA\obik{e‘^*  C el**®-***'*-"*'!.)! 

^  k=l/=l  ® 


The  last  two  terms  of  the  expression  of  Xi  represent  the  crosstalk  in  channel  0.  The 
statistics  of  these  terms  are  extremely  difficult  to  obtain  (if  possible).  Consequently,  for 
tractable  analysis,  Yi  is  modeled  as  a  Gaussian  random  variable.  Such  a  Gaussian 
approximation  has  also  been  used  in  Ref.  8  with  a  different  mathematical  approach. 
Gaussian  approximations  are  commonly  used  to  obtain  the  bit  error  probability  for 
lightwave  systems  when  the  exact  statistics  of  the  decision  variable  cannot  be  analytically 
obtained  [Refs.  13, 17  -18].  From  equation  (6),  the  mean  of  Xi  conditional  on  a  given  bit 
panem  bi  is  given  by 
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(7) 


X*(bi)  = 


From  equation  (7)  and  equation  (AS)  of  Appendix  A,  the  conditional  variance  Ox.  (bi]  of 
Xi  can  be  calculated  as  follows: 


OXi  (bi)  =  JcJ  A Vb?obU2icv/ 4it^(v* + Sj ) 

k>l 

- - ^[2  jcve"^*“(2  jc5k  sin  2Jt5k  -  2icvcos  2jt5k ) 

16it^(v^  +  5k) 

+2it5jte"^’^  (2  JtSk  cos  2ic5^  +  2)tvsin  2jc5k ) + 4ji^  -  4jt^8k]| 


(8) 


where  v  =  pT  and  6^  =  The  parameter  P  is  the  laser  linewidth.  In  summaiy, 

the  Gaussian  approximation  allows  the  decision  variable  Yj  in  equation  (2)  to  be 
considered  as  a  Gaussian  random  variable  with  conditional  mean  Yi(bi)  =  Xi(bi)  in 

equation  (7)  and  conditional  variance  <TYi(^i)  given  by 

<JYi(bi)  =  <yx|(bi)+aw,(bi)+Ow*+ON-  (9) 


For  a  threshold  a,  the  conditional  bit  error  probability  given  bit  patterns  bi  =  (boo. 
bii,-,  biM)  and  b!  =  (bjo.  bn,-,  bjM)  is  [Ref.  19] 


PbC’u.-.biM)  =  |Po(bi')+ jP,(bl) 


(10a) 


where 


(10b) 
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and  eifc(  • )  is  defined  as 


erfc(a)  =  ^  e"**dx .  (lOd) 

The  bit  error  probability  P),  is  obtained  by  taking  the  expectation  of  PbO>ii.*’‘>bjM)  with 
respect  to  the  bit  patterns  (b^i,—,  by^).  Since  there  are  2^  such  patterns,  the  bit  error 
probability  becomes 

where  the  summation  is  over  all  2^  patterns  (l>,i,bi2,‘**.l^y^)-  The  optimal  threshold  that 
minimizes  the  bit  error  probability  is  the  value  that  satisfies  equation  (12). 
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1 

Oy.(b?) 


-{Yi(5l)-«f/2o*.(b|) 


=  0. 


(12) 


Equation  (12)  is  obtained  by  setting  dP(,/da  to  zero. 

In  the  case  when  all  laser  sources  have  the  same  mean  wavelength  but  are 
uncorrelated,  the  above  results  apply  by  setting  5^  =  0  with  k  =  1,  2,—,  M  in  equation 

(8).  Furthermore,  when  all  laser  sources  are  identical  such  that  all  sources  have  the  same 
wavelength  and  phases  noise  process  except  for  the  random  initial  phases,  then  Cx.  (bi) 


in  equation  (8)  reduces  to 


«>i.(Si)  =  iicjAVb?„bL 


(13) 


B.  NUMERICAL  ANALYSIS 

In  this  section  numerical  results  are  presented  for  a  system  with  a  bit  rate  of  500 
Mb/s.  The  responsitivity  of  the  photodiode  is  taken  to  be  0.5  and  the  laser  extinction 
ratio  is  1/20.  The  dark  current  Idk^lO  nA.  The  effective  noise  temperature  for  a  low 
noise  amplifier-integrate  and  dump-slicer  receiver  is  180  K.  Furthermore  assuming  a 
matched  load  of  Rl  =  50  the  post  detection  thermal  noise  current  spectral  density  is 
No  =  2kTo/RL  =  10-22  a2/Hz  [Ref.  15  and  Ref.  20]  where  k  =  1.38  x  10-23  j/k  is  the 
Boltzmann's  constant. 

This  section  is  broken  down  into  the  three  cases  of  laser  sources  and  the  power 
penalty  plot  summaries.  The  data  collected  are  the  results  of  the  MATLAB  model 
(Appendix  C)  of  the  optical  chip  interconnect  for  synchronous  transmission.  The  first 
case  of  laser  source  occurs  when  all  channels  operate  with  independent  laser  sources. 
Case  two  occurs  when  all  laser  sources  have  the  same  mean  wavelength  but  have 
different  noise  processes.  Finally,  in  case  three,  all  laser  sources  are  identical  with  the 
exception  of  the  initial  phase.  These  three  cases  are  summarized  in  the  three  power 
penalty  plot  summaries  for  a  bit  error  probability  Pb  of  10-i3_ 

1.  Case  I 

Figures  5-6  show  the  bit  error  probability  Pb  versus  the  received  peak  power 
A2/2  for  various  levels  crosstalk  from  two  adjacent  channels  relative  to  that  of  a  single 
channel  operation  (zero  crosstalk).  The  laser  linewidth-bit  rate  ratio  and  frequency 
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spacing-bit  rate  ratios  (v,5i,62)  were  taken  to  be  (0.1,0.3,0.3)  and  (0.1.0.7,0.7), 
respectively.  As  seen  in  Fig.  5,  the  bit  error  rate  floor  exists  around  for -20  dB 
crosstalk  irrespective  of  the  received  peak  power.  Qosstalk  levels  must  be  less  than  -26 
dB  for  a  power  penalty  of  1  dB  or  less  at  Pb  = 


Figure  5;  Pb  vs.  A^/2  as  a  function  of  crosstalk  levels  with  (v, 81,62)  =  (0.1, 0.3, 0.3). 


From  Fig.  6,  it  is  seen  that  by  increasing  the  frequency  spacing,  the  crosstalk 
level  can  be  reduced  to  less  than  -23  dB  for  1  dB  or  less  in  power  penalty  at  Pb  =  lO*^^. 
Figures  7-8  show  the  normalized  optimum  threshold  versus  the  received  peak  power  as  a 
function  of  crosstalk  levels.  The  optimum  threshold  decreases  with  increasing  crosstalk. 
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Figure  6:  Pb  vs.  AV2  as  a  function  of  crosstalk  levels  with  (v, 81,52)  =  (0.1, 0.7, 0.7). 


Figure  7:  Normalized  optimal  threshold  a  vs.  as  a  function  of  crosstalk  levels 

with  (v,5i,52)  -  (0.1,0.3,0.3). 
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Figure  8:  Normalized  optimal  threshold  a  vs.  AV2  as  a  function  of  crosstalk  levels 

with  (v^SnSj)  =  (0.1,0.7,0.7). 


Figure  9:  Pb  vs.  as  a  function  of  citsstalk  levels  with  (vtSjfSj)  =  (1,03,0.3). 
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Figure  10:  Pb  vs.  kVl  as  a  function  of  crosstalk  levels  with  (v,5i,52)  =  (5,0J,0J). 


Figure  9  shows  the  results  for  (v,5i,52)  *  (1,0.3, 0.3).  It  is  seen  that  a  laser 
with  a  larger  laser  linewidth-bit  rate  ratio  improves  the  performance.  For  1  dB  or  less  in 
power  penalty  at  Pb  =  lO*'^,  the  permitted  crosstalk  level  is  less  than-23  dB  instead  of 
-26  dB  as  in  Fig.  5.  When  the  normalized  linewidth  is  increased  to  5  as  in  Fig.  10,  there 
is  less  than  1  dB  power  penalty  at  crosstalk  levels  less  than  -17  dB.  This  happens 
because  only  a  portion  of  the  crosstalk  energy  fall  within  the  detection  bandwidth  of  SOO 
MHz.  This  result  encourages  the  use  of  lasers  with  a  large  linewidth  as  long  as  the 
waveguide  bandwidth  is  larger  than  the  signal  spectrum.  When  the  signal  spectrum 
broadened  by  the  laser  phase  noise  appix)aches  the  waveguide  bandwidth,  loss  in  signal 
power  begins  to  occur  and  performance  deteriorates  rapidly. 

For  comparison,  this  model  was  expanded  to  include  the  cases  for  three  and 
four  adjacent  channels  for  synchronous  transmission.  Figures  11  and  12  show  the 
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peiformance  for  the  cases  of  three  and  four  adjacent  channels  in  a  nonplanar  waveguide 
network  with  (v,6i.62,53)=(0.1.0.3.0.3.0.3)  and  (v.5,.82,53,84)=(0.1,0.3,0.3.0.3.0.3). 

respectively.  For  three  adjacent  channels,  a  crosstalk  level  of  less  than  -29  dB  is 
required  for  a  power  penalty  of  1  dB  or  less  at  as  compared  to  a  crosstalk  level 

of  -26  dB  for  the  case  of  two  adjacent  channels  as  shown  in  Fig.  5.  The  bit  error  floor 
for  three  adjacent  channels  exists  at  10*^  for  a  -20  dB  crosstalk  level  irrespective  of  the 
received  peak  power.  In  Fig.  12  for  four  adjacent  channels,  a  crosstalk  level  of  less  than 
approximately  -30  dB  is  required  for  a  power  penalty  of  1  dB  or  less  at  Pb  ®  lO'*^. 
Here,  the  bit  error  floor  exist  around  10*^  for  -20  dB  crosstalk  level.  The  differences 
between  Figs.  11  and  12  are  not  as  dramatic  as  those  seen  between  Figs.  11  and  5  . 
Figures  13  and  14  show  subtle  changes  in  the  normalized  optimal  threshold  in 
comparison  to  Fig.  5.  There  appears  to  be  a  very  slight  change  in  the  normalized  optimal 
threshold  as  the  number  of  adjacent  channels  increases. 


Figure  11:  Pb  vs.  AV2  as  a  function  of  crosstalk  levels  with  (v, 81,82,83)  s 

(0.1,0.3,0.3,0.3) 
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Figure  13:  Normalized  optimal  threshold  a  vs.  kVl  as  a  function  of  crosstalk  levels 

vrith  (v,5j,S2y53)  =  (0.1, 0.3, 0.3, 0.3) 

20 


Figure  14:  Normalized  optimal  threshold  a  vs.  AV2  as  a  function  of  crosstalk  levels 
with  (v, 81,52,53,54)  s  (0.1,03,OJ,03,OJ). 


2.  Case  II 

Figure  IS  shows  the  performance  of  case  two  where  all  laser  sources  have  the 
same  mean  wavelength  but  have  different  noise  processes  for  (v,  61,52)  =  (0.1, 0,0).  A  bit 

error  rate  floor  exists  at  Pb  =  10'^°  for- 20  dB  crosstalk  level.  The  permitted  crosstalk 
level  is  less  than  -29  dB  for  a  power  penalty  of  1  dB  or  less  at  Pb  =  lO"^^.  In  general  the 
performance  of  case  two  is  always  worse  than  that  of  case  one  given  the  same  v  and 
crosstalk  level.  It  is  obvious  that  the  crosstalk  effect  is  reduced  by  using  laser  sources 
with  different  wavelengths. 

3.  Case  m 

Figure  16  shows  the  performance  of  case  three  when  all  laser  sources  are 
identical  with  the  exception  of  the  initial  phases.  The  performance  is  slightly  worse  than 
that  in  Figure  15. 
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Probability  of  Bit  Enor 


Figure  ]5:  Pb  vs.  AV2  as  a  function  of  crosstalk  levels  with  (v,5it52)  s  (0.1, 0,0). 


Figure  16:  Pb  vs.  as  a  function  of  crosstalk  levels  with  all  laser  sources  having 

identical  phase  noise  except  for  initial  phases. 
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4.  Power  Penalty  Plot  Summary 

:  results  for  two  adjacent  channels  with  synchronous  transmission  are 
summarized  in  Figs.  17, 18  and  19  as  power  penalty  relative  to  a  single-channel  operation 
versus  crosstalk  levels,  versus  the  frequency  spacing-bit  rate  ratios  5i  -  Sj,  and  versus 

the  normalized  linewidths,  respectively.  All  results  were  taken  at  Pb  =  lO*^^.  Figure  17 
shows  that  as  the  crosstalk  level  increases,  so  does  the  power  penalty.  It  also  shows  that 
the  smaller  the  normalized  linewidth,  the  more  sensitive  the  system  is  to  smaller  crosstalk 
levels.  Figure  18  shows  that  as  the  power  penalty  increases,  the  normalized  frequency 
spacing  decreases.  Finally,  Fig.  19  suppmts  the  results  in  Fig.  17.  For  normalized 
linewidths  less  than  S,  the  power  penalty  dramatically  increases. 


Figure  17:  Power  penalty  vs.  crosstalk  level  for  a  normalized  frequency  spacing  of 

5|  —  §2  ^  03. 
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Power  Penalty  (dB)  j  n  \  Power  PcraltyCdB) 


18:  Power  penalty  vs.  normalized  frequency  spacing  (5}  =  82)  for  a 
normalized  linewidth  of  v  =  0.1. 


Figure  19:  Power  Penalty  vs.  normalized  linewidth  for  5i  =  62  =  0.3. 
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C.  DISCUSSION 


This  chapter  presented  a  mathematical  ^mework  to  analyze  the  performance  of 
synchronous  optical  chip  interconnects  in  terms  of  the  bit  error  probability  versus  the 
received  power  as  a  function  of  the  crosstalk  level,  frequency  spacing,  and  laser 
linewidth.  This  analysis  can  handle  any  number  of  adjacent  channels.  The  conclusion 
drawn  from  this  investigation  is  that  adjacent  channels  must  use  laser  sources  of  different 
wavelengths  to  reduce  the  effect  of  crosstalk.  Laser  sources  with  a  large  linewidth  also 
help,  as  long  as  the  waveguide  bandwidth  is  much  larger  than  the  resulting  signal 
spectrum.  In  fact,  this  is  the  only  way  to  reduce  the  effect  of  a  given  crosstalk  level  when 
all  laser  sources  have  the  same  mean  wavelength.  When  all  laser  sources  are  locked  to 
the  master  source,  the  performance  depends  explicitly  on  a  crosstalk  level  given  a 
received  peak  power. 

Similar  conclusions  appear  in  Ref.  8  via  a  different  mathematical  approach  that 
applies  to  two  adjacent  channels  only.  The  mathematical  ft-amework  presented  in  this 
chapter  will  be  applied  to  the  asynchronous  case  in  the  next  chapter. 
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ffl.  ASYNCHRONOUS  OPTICAL  CHIP  INTERCONNECTS 


A.  ANALYSIS 

As  in  the  synchronous  transmission  analysis,  for  mathematical  convenience,  the 
complex  envelope  notation  of  a  real  signal  is  adopted.  Thus,  for  a  given  transmitted  bit 
bjo  of  a  given  channel  0  whose  laser  phase  noise  process  is  6o(t),  the  signal  at  the  input 

of  the  photodiode  is  designated  as 

k=l 

k=l 

where  A  is  the  CXDK  signal  amplitude;  bik..i  and  bito  are  the  previous  and  the  present  bit 
in  channel  k  relative  to  bit  bjo;  6^.(t)  is  the  laser  phase  noise  process;  and  ©j^,  ()>^  _|  and 
<{>)(  0  are  the  frequency  spacing  and  the  initial  phase  differences  between  channels  k  and 
0,  respectively.  The  parameter  represents  the  uniformly  distributed  random  delay 
between  bio  and  bik.o.  The  function  u(t)  denotes  the  unit  step  function.  The  parameter 
represents  the  coupling  from  channel  k  to  channel  0.  The  summation  term  in  equation  (1) 
thus  represents  the  crosstalk  from  M  adjacent  channels  into  channel  0.  In  equation  (1),  all 
waveguides  have  the  same  attenuation.  For  the  case  of  M  =  2  and  synchronous  bit 
transmission,  that  is,  only  the  two  nearest  adjacent  channels  are  considered,  equation  (1) 
reduces  to  the  result  in  equation  (2.4)  in  Ref.  8.  In  this  case  Cj.  =  b/‘v1-2B^  ,  k=l,2 
where  B  is  the  total  power  coupled  from  a  waveguide  into  its  adjacent  waveguide.  In 
practice,  bjjc,.!  and  bac^o  *  0,1,— ,M  are  not  necessarily  equal  to  zero  or  one.  Let  k  = 
0,1,'**,M  be  the  extinction  ratio  of  the  laser  of  channel  k  (defined  as  the  ratio  of  the 
transmitted  power  of  the  logical  zero  to  that  of  a  logical  one).  Then 
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bok.o.bok.-i  *  ‘>j^k/{^'*‘h)  logical  zero  and  bi^o«bik,_i  *  logical 

one.  Hereafter,  ri^'s  are  assumed  identical  for  ail  channels. 

Let  Rsn^q/hf  be  the  photodiode  responsitivity  [Ref.  14  -  Ref.  15]  where  Oe  S  1  is  the 
quantum  efficiency,  q  is  the  electron  charge  (1.6x10-^^  C),  h  is  Planck's  constant  (6.626  x 
10-34  f  ij.  frequency.  The  output  of  the  photodiode  is 

R|si(t)j^  +  Wi(0+ '^dk(0  where  w,(t)  is  the  shot  noise  generated  by  the  photodiode  and 
w^(t)  is  the  dark  current  noise.  The  output  of  the  photodiode  plus  the  post-detection 
thermal  noise  n(t)  is  integrated  by  the  integrare-and  dump  ffiter  with  the  normalization 
constant  R  resulting  in  the  decision  variable  Y|  as  follows: 


=  Xi-^W,  +  Wdlt  +  N 


(15a) 


where 


(15b) 

W,=i)„w.(.)d. 

(15c) 

W*=ijJwik(t)dt 

(15d) 

(15e) 

Since  n(t)  is  a  zero  mean  Gaussian  process  with  spectral  density  No,  the  Gaussian  random 
variable  N  also  has  zero  mean  and  its  variance  is  given  by 


(16) 
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On  the  other  hand,  the  shot  noise  Wg(t)  is  a  non-stadonaiy  process  since  the  envelope  of 
the  signal  at  the  input  of  the  photodiode,  namely  isi(t)l,  is  tinK-dependent.  Therefore,  the 
shot  noise  Wg(t)  can  be  nxxieled  as  a  zero  mean  wide-sense  stationary  Gaussian  process 
whose  spectral  density  Wo(bi),  given  a  bit  pattern  bi  « (bio,bii.-i.  •.biM.-i.bii.o,  •,biM,o).  is 

proportional  to  the  conditional  mean  of  the  squared  envelope  of  the  input  signal.  In  other 
words,  E||si(t){^|bi|  approximates  lsi(t)l^  over  a  bit  time  T.  Based  on  this  approximation, 

the  shot  noise  spectral  density  can  be  obtained  as  follows  [Ref.  17] 

W„(bi)  =  qRE{|Si(lf|bi} 


From  equation  (ISc)  and  equation  (17),  the  conditional  variance  of  Wg  given  a  bit  panem 
bi  is 


(18) 


The  dark  current  noise  spectral  density  function  is  ql^^  where  Idk  is  the  dark  current.  The 
variance  of  the  dark  current  noise  is  =  Tqldk/R^-  The  random  variable  Xi  in 

equation  (ISb)  consists  of  the  signal  term,  the  signal-crosstalk  terms,  and  the  crosstalk- 
crosstalk  terms. 

The  statistics  of  the  crosstalk  terms  are  extremely  difficult  to  obtain  (if  possible). 
Consequendy,  for  tractable  analysis,  Yj  is  modeled  as  a  Gaussian  random  variable.  Such 
a  Gaussian  approximation  has  also  been  used  in  Ref.  8  for  synchronous  transmission  with 
a  different  mathematical  approach.  Gaussian  approximations  are  commonly  used  to 
obtain  the  bit  error  probability  for  lightwave  systems  when  the  exact  statistics  of  the 
decision  variable  cannot  be  analytically  obtained  [Refs.  13,17  -18].  From  equation  (ISb), 
the  mean  of  Xi  conditional  on  a  given  bit  pattern  bi  is  given  by 
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(19) 


X,(bi)  =  E{Xi|bi} 

-  |a*TT)5,  +iXciAhl>L..,  +^XCkA^i.o- 

From  equation  (18)  and  equation  (B5)  of  ^pendix  B,  the  conditional  variance  (bi) 

of  Xi  can  be  calculated  as  follows: 


®Xi  (bi)  =  X^kA‘^^bfo(bi..i  +  b5^o){xv/ 4jc^(v^ + ) 

ksl 

64ji®(v2  +  85) 

_g-2itv  ^2jcvsin  2jc8k  +  2118^  cos2ii8^ )] 

44rc^(Sk  -y^)[2jcv-4)c^(v^  +  8j) 

^-2icv(2jjg^  sin  2jt8k  -  2nvcos2ic8^  )]| 


where  v  =  PT  and  8^  -  oi)^T/27c  .  Here  3  is  the  laser  linewidth.  In  summary,  the 

Gaussian  approximation  allows  the  decision  variable  Yj  in  equation  (2)  to  be  considered 
as  a  Gaussian  random  variable  with  conditional  mean  Yi(bi)  >  Xi(bi)  in  equation  (19) 
and  conditional  variance  Oy;  (bi)  given  by 

OYi(bi)*<yxj(bi)+o^^(bi)+aw^ +On-  (21) 

For  a  threshold  a,  the  conditional  bit  error  probability  given  bit  patterns  bf  =  (boo. 
bii.-i.'".  biM,-i.  bii.0,— .biM.o)  and  bj  =  (bio,  bii,.i,— ,  biM,-i,  bii,o,— ,biM.o)  is  [Ref.  20] 

Pb(bii -i»—»biM -i.bii^O’***’biM,o)  “  2 2 (^^) 

where 
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where 


and  erfc(  * )  is  again  defined  (as  in  equation  (lOd))  as 


(22b) 


(22c) 


(22d) 


The  bit  error  probability  Pj,  is  obtained  by  taking  the  expectation  of  Pb(bii,.i,— ,  biM.-i. 
bii.o.-.biM.o)  with  respect  to  the  bit  patterns  (bii..i,—,  bjM.-i.  bii,o,-,biM.o)*  Since  there  are 
22M  5uch  patterns,  the  bit  error  probability  becomes 


*••••  tiikl.o) 


^b{bii,-i  »•••»  bjM.-i  >  bji,o*““>  bjM.o )  (23) 


where  the  summation  is  over  all  2^^  patterns  (bii..i,—,  biM.-i,  bii,o,— ,biM,o).  The  optimal 
threshold  that  minimizes  the  bit  error  probability  is  the  value  that  satisfies  equation  (24). 

G4) 


aY,(b?) 


OYj(b!) 


j  ,-[Y.(a)-.]’/244fi) 


=  0 


Equation  (24)  is  obtained  by  setting  dP^/da  to  zero. 

In  the  case  when  all  laser  sources  have  the  same  mean  wavelength  but  are 
uncorrelated,  the  above  results  apply  by  setting  =  0  ,  k  =  1,  2,  -.M  in  equation  (21). 
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Funhennore,  when  all  laser  sources  are  identical  such  that  all  sources  have  the  same 
wavelength  and  phases  noise  process  except  for  the  random  initial  phases,  then  Ox,(^i) 

in  equation  (20)  reduces  to 

°x,(bi)  =  1^  iciA*T»bf„(bi.,  +  bi.,,).  (25) 


B.  NUMERICAL  ANALYSIS 

In  this  section,  numerical  results  are  presented  for  a  system  with  a  bit  rate  of  S(X) 
Mb/s.  The  responsitivity  of  the  photodiode  is  taken  to  be  0.5  and  the  laser  extinction 
ratio  is  1/20.  The  dark  current  =  10  nA.  The  effective  noise  temperature  of  a  low 
noise  amplfier-integrate  and  dump-slicer  receiver  is  180  K.  Furthermore,  assuming  a 
matched  load  of  =  50  fl,  the  post-detection  thermal  noise  current  spectral  density  is 
No  =  2kTo/RL  =  10-22  a2/Hz  [Ref.  15  and  Ref.  20],  where  k  =  1.38  x  10*23  J/K  is  the 
Boltzmann's  constant. 

This  section  is  broken  down  into  the  three  cases  for  the  probability  of  bit  error  plots 
and  the  power  penalty  plot  summaries.  The  power  penalty  plots  compare  the  case  of  two 
adjacent  channels  having  synchronous  transmission  with  the  case  of  two  adjacent 
channels  having  asynchronous  transmissicn.  The  cases  summarized  in  the  three  power 
penalty  plots  are  for  a  bit  error  probability  Pb  of  10*^  5.  The  data  collected  is  the  result  of 
the  MATLAB  model  of  the  optical  chip  interconnect  for  asynchronous  transmission 
(Appendix  D)  and  synchronous  transmission  (Appendix  C).  The  first  case  of  laser  source 
occurs  when  all  channels  operate  with  independent  laser  sources.  Case  two  occurs  when 
all  laser  sources  have  the  same  mean  wavelength  but  have  different  noise  processes,  and 
case  three  occurs  when  all  laser  sources  are  identical  with  the  exception  of  tlie  initial 
phases. 
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1.  Case  I 

Figures  20-21  show  the  bit  -  Probability  P|,  versus  the  received  peak  power 
A2/2  for  various  levels  of  crosstalk  from  two  adjacent  channels  relative  to  that  of  a  single 
channel  operation  (zero  crosstalk)  and  with  laser  linewidth-bit  rate  ratio  and  frequency 
spacing-bit  rate  ratios  (v,5i,52)  taken  to  be  (0.1,0.3,0.3)  and  (0.1, 0.7 ,0.7),  respectively. 

From  Fig.  20,  the  bit  error  rate  floor  exists  around  Pj,  =  lO'*^  for -20  dB  crosstalk 
irrespective  of  the  received  peak  power.  Crosstalk  levels  must  be  less  than  -26  dB  for  a 
power  penalty  of  1  dB  or  less  for  Pb  =  10**5.  From  Fig.  21,  it  is  seen  that  by  increasing 
the  frequency  the  frequency  spacing  the  crosstalk  can  be  reduced  to  less  than  -23  dB  for 
1  dB  or  less  in  power  penalty  at  Pj,  =  10**5. 


Figure  20:  Pb  vs.  AV2  as  a  function  of  crosstalk  levels  with  (v,5|,52)  =  (0.1,0.3,0.3)> 
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Figure  21:  Pb  vs.  A^/l  as  a  function  of  crosstalk  levels  mth  (v,5i,S3)  =  (0.1, 0.7, 0.7). 


Figure  22:  Normalized  optimal  threshold  a  vs.  AVZ  as  a  function  of  crosstalk  levels 

with  (v, 81,52)  =  (0.1, 0.3, 0«3). 
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Figure  23:  Normalized  optimal  threshold  a  vs.  AV2  as  a  function  of  crosstalk  levels 

with  (v,5„52)  a  (0.1, 0.7, 0.7). 


Figures  22-23  show  the  normalized  optimum  threshold  versus  the  received 
peak  power  as  function  of  crosstalk  levels.  The  optimum  threshold  decreases  with 
increasing  crosstalk.  Figure  24  shows  the  results  for  (v, 81,52)  =  (1,0.3, 0.3).  It  is  seen 

that  a  laser  with  a  larger  laser  linewidth-bit  rate  ratio  improves  the  perfoimance.  For  0.7 
dB  or  less  in  power  penalty  when  Pj,  =  10*1^,  the  permitted  crosstalk  level  is  less  than 
—23  dB  instead  of  -26  dB  as  in  Fig.  20.  When  the  normalized  linewidth  v  is  increased 
to  5  as  in  Fig.  25,  there  is  less  than  0.7  dB  power  penalty  for  crosstalk  levels  less  then 
- 17  dB.  This  happens  because  only  a  portion  of  the  crosstalk  energy  falls  within  the 
detection  bandwidth  of  500  MHz.  This  result  encourages  the  use  of  lasers  with  a  large  v 
as  long  as  the  waveguide  bandwidth  is  larger  than  the  signal  spectrum.  When  the  signal 
spectrum  broadened  by  the  laser  phase  noise  approaches  the  waveguide  bandwidth,  loss 
in  signal  begins  to  occur  and  performance  deteriorates  rapidly. 
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Figure  24:  Pb  vs.  kVl  as  a  function  of  crosstalk  levels  vrith  (v,5|,52)  =  (1,03,03). 


Figure  25:  Pb  vs.  A2/2  as  a  function  of  crc»stalk  levels  vrith  (v, 5], 82)2  (5,03,03). 


2.  Case  II 

Figure  26  shows  the  performance  of  case  two  when  all  laser  sources  have  the 
same  mean  wavelength  but  with  different  noise  process  for  (v.5i,52)  » (0.1. 0, 0).  A  bit 

error  rate  floor  exists  at  P),  =  lO-^^  irrespective  of  the  received  peak  power.  The  permitted 
crosstalk  level  is  less  than  -26  dB  for  a  power  penalty  of  1  dB  or  less  at  P), «  lO*^^.  In 
general,  the  performance  of  case  two  is  always  worse  than  that  of  case  one  given  the 
same  and  crosstalk  level.  It  is  obvious  that  the  crosstalk  effect  is  reduced  by  using  laser 
sources  with  different  wavelengths. 


Figure  26:  Pb  vs.  AV2  as  a  function  of  crosstalk  levels  with  (v.SiySj)  s  (0.1,0,0). 
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3.  Case  III 


Figure  27  shows  the  performance  of  case  three  when  all  lasers  are  identical 
with  the  exception  of  the  initial  phases.  The  performance  is  slightly  bener  than  that  in 
Fig.  26.  These  results  relate  directly  back  to  the  derivation  of  equation  (25)  and  are 
opposite  to  those  in  Case  III  for  synchronous  transmission. 


-60  -38  -56  -54  -52  -30  -48  -46  -44 

Peak  Powei<dBW) 

Figure  27:  Pb  vs.  AVI  as  a  function  of  crosstalk  levels  with  all  laser  sources  having 
identical  phase  noise  processes  except  for  the  initial  phases. 

4.  Power  Penalty  Plot  Summary 

Figure  28  shows  the  power  penalty  versus  crosstalk  level  for  two  interfering 
channels  with  a  normalized  frequency  spacing  of  0.3  as  a  function  of  the  normalized  laser 
linewidth.  The  power  penalty  for  the  synchronous  transmission  where  5]^  &  0.3  in 
equation  (1)  is  also  plotted  for  comparison.  In  general,  the  power  penalty  is  smaller  for 
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asynchronous  transmission  when  the  linewidth  is  less  than  3  for  crosstalk  levels  up  to 
-20  dB.  For  larger  laser  linewidths,  the  power  penalty  for  synchronous  transmission  is 
slightly  less.  Figure  29  shows  the  power  penalty  versus  normalized  frequency  spacing 
for  two  interfering  channels  with  a  normalized  linewidth  of  0.1  as  a  function  of  the 
crosstalk  level.  In  general,  the  power  penalty  for  asynchronous  transmission  is  less  than 
that  of  synchronous  transmission.  Hgure  30  shows  the  power  penalty  for  two  interfering 
channels  versus  the  normalized  laser  linewidth  with  a  normalized  frequency  pacing  of 
0.3,  as  a  function  of  the  crosstalk  level.  In  general,  the  power  penalty  for  asynchronous 
transmission  for  normalized  laser  linewidths  less  than  3  is  less  than  that  for  synchronous 
transmission. 


Figure  28:  Power  Penalty  vs.  crosstalk  level  for  a  normalized  frequency  spacing  of 
5}  =  82  =  03  for  asynchronous  and  synchronous  transmission. 
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Figure  29:  Power  penalty  v&  normalized  frequency  spacing  (5|  -  62)  for  a 
normalized  linewidth  of  v  =  0.1  for  asynchronous  and  synchronous  transmission. 


Figure  30:  Power  penalty  vs.  normalized  linewidth  for  81 -bj-  0.3 
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C.  DISCUSSION 

To  conclude,  this  chapter  presented  a  mathematical  framework  to  analyze  the 
performance  of  asynchronous  optical  chip  interconnects  in  terms  of  the  bit  error 
probability  versus  the  received  peak  power  as  a  function  of  the  crosstalk  level,  frequency 
spacing,  and  laser  linewidth.  Furthermore,  the  analysis  can  handle  any  number  of 
adjacent  channels.  The  conclusion  drawn  from  this  investigation  is  that  adjacent  channels 
must  use  laser  sources  of  different  wavelengths  to  reduce  the  effect  of  crosstalk.  Laser 
sources  with  a  large  linewidth  also  help,  as  long  as  the  waveguide  bandwidth  is  much 
larger  than  the  resulting  signal  spectrum.  In  fact,  this  is  the  only  way  to  reduce  the  effect 
of  a  given  crosstalk  level  when  all  laser  sources  are  locked  to  a  master  source.  The 
performance  depends  explicitly  on  the  crosstalk  level  given  a  received  peak  power. 
Finally,  as  seen  by  Figs.  28, 29.  and  30,  asynchronous  optical  chip  interconnects  perform 
better  than  their  synchronous  counterparts. 
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VL  CONCLUSION 


This  thesis  derived  a  crosstalk  nxxlel  for  an  optical  chip  interconnect  that  determines 
the  crosstalk-induced  system  penalty  caused  by  the  coupling  between  two  or  more 
adjacent  channels  in  a  single-mode  waveguide  array.  MATLAB  version  3.Sk  was  useful 
in  the  development  and  the  verification  of  the  nrodel.  The  resulting  perfmmance  curves 
and  power  penalty  plot  summaries  demonstrate  the  effectiveness  of  the  nxxlel  for  both 
synchronous  and  asynchronous  transmission.  In  general,  the  asynchronous  case  proved 
to  be  more  sensitive  than  the  synchronous  case  fen  crosstalk  levels  less  than  -20  dB  and 
normalized  laser  linewidths  less  than  2.  These  results  follow  those  in  Ref.  9  where  it  was 
determined  that  the  sensitivity  of  an  asynchronous  receiver  is  only  about  0.S  dB  less  than 
its  synchronous  counterpart 

Study  is  currently  underway  to  determine  the  effect  on  this  crosstalk  model  of 
shutting  down  one  or  more  of  the  adjacent  channels  during  transmission.  The  results 
from  this  study  should  further  verify  the  effectiveness  of  this  model. 
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APPENDIX  A  -  DERIVATION  OF  THE  CONDITIONAL 
MEAN  SQUARE  VALUE  OF  X,  GIVEN  bi  FOR 
SYNCHRONOUS  TRANSMISSION 


The  conditional  mean  square  value  of  Xj  given  bi  is  important  in  the  derivation  of  the 
conditicmal  variance  (bi)  of  Xj  (equation  (8)).  The  derivation  of  the  conditicxnal  mean 

square  begins  by  finding  the  expectation  of  Xf  (sec  equatiem  (6))  given  bi. 


E{x?|bi}  = 


a  E 


+  2  X^kA^biobacI^  Je  dt 

+1 X  X  CkQA^bikbio  fe  -«>i)i+0k  -♦ildtj 

^  k=ir=i  ® 

+|XCiA''Tbfobn[/J' 
+“XXCkC/A'*TbJobikbi,rcjl®k(‘>-®'(‘H®k-‘0/)t+«k-Oi]dt 

2  k*l/=l  *'® 

+71  ickC<A'bfob»b  J r  (•)-».-♦»  I* 

^  k*l/=l  L  ® 

+ilSic.CjC,A‘bi„bi,bnbJjV»«<‘)-».(')-».-».ld, 

■^n=lk=l/*l  ‘-•'W 
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T  -jteo(t)-e.(t)-«).t-4.]jjjj|'*’gXek(T)-e^(t)+K-®/h+^k-0/]dx 


+L  c 


1  M  M  M  M 

+7 1 1  IIC„C„CkC,ASi„bi„bikbu 


mslnslksl/sl 


xj; 


.Tcjl®k  )*+♦»  -♦■Idt 


iA'rVo+iX'^A'Vbfobi 

+TScjA*b?ob5, 

4  j 

■fa  g- jo)i  (t-t)£|g  j(eo(t)-ek  (t)-eo(t)+ek  (t))|jj 
<Jo  g  j«>k  (‘“'*)E|g*  jl®o(l)-®k  (0-®o('*)+®k  (■*)]  jdtdx  j 


■^TiicicjAVbib^^ 
‘♦k=ls 


(A.1) 


In  the  derivation  of  equation  (A.l),  all  the  initial  phase  differences,  <j>if,  k  =  1, 

M  are  uniform  variables  over  (0,  2k). 

Let  v  =  0o(t)-0k(t)-[0o('c)-®kWl-Thc^®®®’’P^“®"o^®®  ®k(0*  k  =  0,  1,—,  M 
is  characterized  by  a  Wiener  process  [Ref.  17  -  Ref.  18]  such  that  d0k(t)/dt  =  2Kp(t) 
where  p(t)  is  a  zero  mean  white  Gaussian  process  of  PSD  P/2k  Hz  where  P  is  the  laser 
linewidth  (assumed  to  be  the  same  for  all  laser  sources).  The  variance  of  0k(t)  is  2KPt. 
The  process  y  is  a  zero  mean  Gaussian  process  and  can  be  expressed  as 

V  =  2K0Mti)-Hk(ti)]dti  “2K0Ho(ti)-lik(ti)}iti 

(A.2) 
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Therefore,  the  variance  is  given  by 


=  4jt^rf  •^(*i”*2)dtidt2 

=  4nPjt-tl. 

Using  the  fact  that  y  and  *-y  are  both  Gaussian  random  variable  with  zero  mean  and 
variance  47[p|t  -  t|,  the  expectation  value  is  [Ref.  13] 

E{ei»}  =  E{e-i»}  =  e-'iA  =  (A.4) 

Substituting  equation  (A.4)  into  equation  (A.1),  the  conditional  mean  variance  is  simplified 
to 

E{x?|bi}  =  ■jAVbfo  +  iXclAVbfobi^ 

+^S^A'bfob§,  r  cos(i),(t  -  tjdtdi 

^  k=l/=l 

=  iA*T2bf„+||;(^AVb5,bi 

^  '^k=l 

+7  SCkA^b?oba/^(T  -  |u|)c"^"P*“l  cosm^udu 

*^k=l 

^  k=ir*i 

=  iAVb?o+ifc^A'T*b?„bi 

+iiic5ciAWibi 

^k=l/=l 
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+XQAVbfobi 

ksl 


2jcv 


l4it^(v^+5k)  16Ji^(v2+6jf 

>^2jcvc"^*^(2jc5i;  sin2icS)^  -  litvcoslnb^) 


(A.5) 


+2n&^c~^^{2nS^co&2nb^  +2jtvsin2jiOk) 


+4jt 


^^-471^8^]} 


where  V  =  px  and  5^  =  <0|jT/2jc. 
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APPENDIX  B  -  DERIVATION  OF  THE  CONDITIONAL  MEAN 
SQUARE  VALUE  OF  X,  GIVEN  bi  FOR  ASYNCHRONOUS 

TRANSMISSION 


The  conditional  mean  square  value  of  Xi  given  bi  is  important  in  the  derivation  of 
the  conditional  variance  <Tx^(bi)  of  (equation  (20)).  The  random  variable  X^  is 

derived  from  equations  (14)  and  (ISb).  By  taking  the  expectation  of  xf  given  bi  and  x 
where  t  =  ('Ci,X2,...,Xm)  .  the  conditional  mean  square  is 

x[u(t)-u(t-Xij)]dt 

_  u(t  _  )]dt  j 

k=l  ‘- 

x[u(t  - )  -  u(t  “  T)]dt 

+ )  -  u(t  -  T)]dtJ 

M  M 
k=l/=l 

X  ejl®k  (i)-®/(‘)+(®k  -®/)‘+^ik,-i  -♦i/.-i  ] 

Jo 

x[u(t)  -  u(t  -  Tk )  Ju(t)  -  u(t  - 1/  )]dt 

M  M 

+X  X^kC/bik._ibi/,o 

k=ir=i 

X  r  e  j(®‘  ‘)+(®k  -®/  )*+^ik.-i  -♦u.o  1 

Jo 

x[u(t)  -  u(t  -  Xk  )|u(t  -  X/ )  -  u(t  -  T)]dt 
M  M 

+XX^kCrbik_ibi^,o 
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^  f  g- j[ek  ( t)-e<  (i)+(®k  -»/  )»+♦*  .-1  -♦i/.o  ] 

Jo 

x[u(t)  -  u(t  -  tk ) Ju(t  -  ti )  -  u(t  -  T)]dt 
M  M 

+ X  S  ^kQ^ik.O*^i^.O 

ksl/sl 

X  C  -“i  )‘+^ik.o-«  i/  .0 ) 

2 

x[u(t  -  tk )  -  u(t  -  T)Ju(t  -  T^)  -  u(t  -  T)]dt}  . 


(B.l) 


Expanding  the  terms,  the  conditional  mean  becomes 


Efxflbi.ij  =  ■^EjT2bfo  +  XcJbfobi..,Tti+£cJb?obioT(T-ti) 


M 


+Ic5b?ob2., 

k=l 


M 


.2u2  k2 

ik.O 


■XcSbfob 


k=l 


J^e"-'®*  (‘-‘'M®o(*>‘®k(‘>"®‘>^^)'"®k  (‘'lldtdx 

^Wxj 


+  f  J  e 

Jtk  Jtk 
M 


M  M 


+XcJb?obL..iTXk  +  XZcjQbi,.ibJ._iXkX^ 

k=i  k=ir=i 
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k>U-l 

+ ^  Ckbfob^,oT(T  "  “^k ) 

k>l 

+SXc5cJbi..ibJ.otk(T-T,) 

k>l/sl 

+iicic5bi.ob5.o(T-XkXT-t,) 

k*U»l 


Evaluating  the  expectation  of  each  term,  equation  (B.2)  becomes 


E{x?lbi.t}  =  -^{T^bfo  +  lXc^bJjbi.iTtk 


k^\ 

M 

+2XCkb5,bi,oT(T-t,) 

ksl 

+22^]^CkC/b^  _ib5,oXk(T”  Xjc) 

ksl/sl 

+ZZcScJbi,.ibJ._itkT/ 

k®l/=l 

+Zicic5biob5,o(T-  tkXT- 

k=l/=l 

+^Cijbjobik.-i 

k=l 


p  ^jo.k(t-t)£|g-ieo(t)-ek(i)-eo(T)+OkW]|<jtdtj 


M 


+2Ckbioblo 

k-l 


(B.2) 
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(B.3) 


The  conditional  variance  of  Xj  given  by  hi  and  t  is  given  by 


<^,(bi,i)  =  E{x?lbi.i}-[Xi(bi.i)f 

a4  f  M 

=  "j"|SCkbfob^_i 

■^rr  gj«#k(t-t)£j^-Xeo(t)-ek(i)-Oo(t)+ek(‘')]|(jtdtj 

+5^Ckbiob^o 

lc»l 

iCC 

<x  g  j®k  (*“‘*)E|c”j^®®  ^*^0  (■*h®  k  Wljdtdt j| . 


(B.4) 


It  should  be  noted  that  the  derivation  of  equation  (B.4)  depends  on  the  fact  that  all  initial 
phase  differences  ^  with  k  s  1, 2,...,  M  are  uniform  random  variables  over  (0,27c). 

Let  t|r  s  eo(t)  -  0k(t)  -  [OoCt)  -  0k (t)].  The  laser  phase  noise  0k(t)  with  k  *  0, 1,-, 
M  is  characterized  by  a  Wiener  process  [Ref.  17  -  Ref.  18]  such  that  d0k(t)/dt  s  2iqi(t) 
where  p.(t)  is  a  zero  mean  white  Gaussian  process  of  PSD  P/2ic  Hz  where  P  is  the  laser 
linewidth  (assumed  to  be  the  same  for  all  laser  sources).  The  variance  of  0^(t)  is  In^t. 
The  process  v  is  &  zero  mean  Gaussian  process  and  can  be  expressed  as 
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V  *  2icJo^o(ti)-  ^k(t|)K  “  2R0^io(ll)-  ^lk(tl)]dtl 

*  2Ji||[Ho(ti)- ^k(ti)ldti .  (B.5) 

Therefne,  the  variance  is  given  by 

=  “■  ‘'I "  1“!)^“)^“ 

=  4jlP|t-x|. 

Using  the  fact  that  y  and  -y  are  both  Gaussian  random  variables  with  zero  mean  and 
variance  4nP(t-  tl,  the  expected  value  is  [Ref.  13] 

E{e>»}  =  E{e->»}  =  (B.7) 

Substituting  equation  (B.7)  into  (B.4),  the  conditional  variance  obtained  is 

^  lk=l 

+X^^fobL.o  cos©k(t  -  T)dtdx| 

=  (Xk-H)c'^’^'"’cos©kUdu 

^  lk=l  "  * 

Finally,  the  conditional  variance  of  Xj  given  bi  is  obtained  by  averaging  Ox.  (bi.t) 
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over  t  giving  the  end  result  as 


=  XcjA'T’bJ, 

ksl 


1 

64ii*(v^+8j)^ 

x(8it*v5k[2ii5k  -e‘**'(2jcvsin2ii5k  +  27c5k  cos2ji8k)] 

+4jc*(5j-v^) 

x|2icv  -  4k^(v*  +  8k) 

x(2ji5ij  sin27t8ij  -27tvcos2ji8^)])j 


(B.9) 


where  v  =  PT  and  8ij  =  a)yT/2n. 
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APPENDIX  C  -  MATLAB  MODEL  OF  OPTICAL  CHIP 
INTERCONNECT  FOR  SYNCHRONOUS  TRANSMISSION 


Ibis  model  is  for  the  case  of  2  adjacent  channels  with  (v,5|,52)  -  (0.1,0.3,0.3)  for 
synchronous  transmission.  The  MATLAB  version  used  was  version  3.5k.  The  program 
itself  is  a  translation  of  equations  (1)  through  (13)  into  MATLAB  code. 

%  Synchronous  Case:  No=:lE-22 
%  Case  1.3:  nu=.l  and  delta  function  subscript  k  (delk)=.3 
%  th2chl3x.m 
%  19  July  93 
%  Known  constants 

R=.5; 

T=1/(500E6); 

N0=lE-22; 

q*1.6E-19; 

IdkslOE-9;  %  dark  current 

%  Coupling  constants  for  Ck^2 
ckl=0; 

ck2=3.16E-4;  %-35dB 
ck3=6.32E-4;  %-32dB 
ck4=1.26E-3;  %-29dB 
ck5=2.52E-3;  %-26dB 
ck5a«5.011E-3;  %-23dB 
ck6=.01;  %-20dB 

ck7=1.995E-2;  %-17  dB 
ck8=3.981E-2;  %-14  dB 
ck9=7.943E-2;  %-ll  dB 
ckl0=.1585;  %-8dB 

ck=[0  3.16E-4  6.32E-4  l,26E-3  2.52E-3  5.012E-3 .01  1.995E-2  3.981E-2  7.943E-2 
.1585]; 

%  Pattern  bit  values 

%bl=[.2182 .9759 .9759 .2182];  %for  bll.bOl  for  patterns  A3.CJ) 

%b2=[.9759 .2182 .9759 .2182];  %for  bl2,b02  for  patterns  A3.CJ) 

blA=.2182; 

blB=.9759; 

blO.9759; 

blD=.2182; 

b2A=.9759; 

b2B».2182; 

b2C=.9759; 

b2D=.2182; 

b3A=bl/  "2+b2A''2; 

b3Bsblb  2+b2B^2; 

b3C*blC''2+b2C''2; 

b3D=bliy'2+b2D^2; 
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500^^182;  %  approximate  zero 
bl0B.S>759;  %  approximate  1 

%  Determine  k  constant  for  the  variance  of  X 

nu«.l; 

delk».3; 

ka«(2*pi*nu)/(4*pi'^2*(nu'^2+delk'^2)); 

kb=l/(16*piM*(nu''2+delk''2)^2); 

kc*(2*pi*nu*exp(-2*pi*nu))*(2*pi*delk*sin(2*pi*delk)-2*pi*nu*cos(2*pi*delk)); 

kd*(2*pi*delk*exp(-2*pi*nu))*(2'^i*(delk*cos(2*pi*delk)+nu*sin(2*pi*delk))); 

ke=4*pi'^2*(nu'^2-delk'^2); 

k»ka-kb*(kc+kd+ke); 

n«l; 

j»i; 

m*l; 


%  -60dB<*Peak  Power<=-44dB 
%  Peak  power  is  A'^2/2  »  Ps 
for  ss=-60:.2;-44; 
ps(m>=l(yHss/10); 
m=m+l; 
end 

%  The  goal  here  is  to  solve  the  equation  fen-  the  optimal  threshold 
%  vs  p^  power 

alpha=zeros(length(ps),length(ck)); 

PBEszeros(length(ps),len^(ck)); 

for  js  1  :length(ck);  %  coupling  values  loop 
n=l; 

for  n-l  :length(ps);  %  peak  ppower  values 
Ps=ps(n); 

%  Constants 

qcnst=q/(Ps*T*R); 

vnoise=N0/(R''2*Ps^2*T); 

vdark=(q*Idk)/(R^2*Ps''2*T); 

%  Determine  the  sigma  for  each  of  the  0  bit  patterns 

sOA*ssqrt(4*ck(j)*bO(y'2*k*b3A+qcnst*(b(X>'2+ck(j)*b3A)-fvnoise+vdark); 

sOB=s^(4*ckQ)*bO(y'2''‘k*b3B+qcnst*(bO(y^2+ck(j)*b3B)+vnoise+vdark); 

s0C*sqrt(4*ck(j)*b00^2*k*b3C+qcnst*(b0CK^2+ck^)*b3C)+vnoise+vdark); 

s0D=sqrt(4*ck^)*b00^2*k*b3D4qcnst*(b0(y'2+ck5)*b3D)+vnoise+vdaik); 

%  Detmnine  the  variance  for  each  of  the  0  bit  patterns 


sOA2*sOA^2; 

sOB2=sOB'^2; 

sOC2=sOC''2; 

sOD2=sOD^2; 

%  Determine  the  sigma  for  each  of  the  1  bit  patterns 
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slA»sqrt(4*ck(i)*bl(K^2*k*b3A+qcnst*(blCy'2+ck(j)*b3A)+vnoise+vdark); 

slB*s^(4*ck(j)*bl0''2*k*b3B-Hicnst*(bl0^2+ck(j)*b3B)+vnoisc+vdaric); 

slC«s^(4*ck(j)*bl0^2*k*b3C4qcnst*(bl(y^2+ckO)*b3C)+vnoisc+vdaric); 

slD«s^4*ck^)*bl(y'2*k*b3D-Kicnst*(bl(y'2-«:kO)*b3D)+vnoise+vdark); 

%  Determine  the  variance  for  each  of  the  1  bit  patterns 

slA2»slA'^2; 

slB2«slBA2; 

slC2«slC''2; 

slD2»slDA2; 

%  Determine  the  mean  for  each  of  the  0  bit  patterns 

mOA«bO(y'2+ck(j)*b3A; 

mOB=b00^2+ck(j)*b3B; 

m0C»b00^2+ck(j)*b3C; 

m0D=b00^2+ck(j)*b3D; 

%  Determine  the  mean  for  each  of  the  1  bit  patterns 

mlA=bl0^2+ck(j)*b3A; 

mlB=blO^2+ck0)*b3B; 

mlC=blO^2+ck0)*b3C; 

mlD=bl0^2+ck(j)*b3D; 

%Solve  the  equation  F(a)=0  to  optimize  a,  the  threshold.  F(a)sO  comes 
%  from  (dP/(la)»0 

a=0:.001:l' 

FA=(l/s6A)*exp(-((a-mOA).'^2)./(2*sOA2))-(l/slA)*exp(-((mlA-a).'^2)/(2*slA2)); 

FB=(l/s0B)*exp(-((a-m0B).A2)./(2*s0B2))-(l/slB)*exp(-((mlB-a).'^2)/(2*slB2)); 

FC=(l/s0C)*exp(-((a-m0C).''2)./(2*s0C2))-(l/slC)*exp(-((mlC-a).''2)/(2*slC2)); 

FD=(l/s0D)*exp(-((a-m0D).'^2)./(2*s0D2))-(l/slD)*exp(-((mlD-a).''2)/(2*slD2)); 

F=FA+FB+FC+FD; 

[YJ]=min(abs(F)); 

alpha(nj)=a(l); 

aasalpha(nj); 

%  Determine  the  probability  of  bit  error  for  each  optimal  threshold 
%  ovCT  all  bit  patterns  A.P  CD 
pOA(n)=aerf((alpha(nj)-mOA)/sOA4nf); 
pOB(n)serf((alpha(nj)-mOB)/sOB4nf); 

pOC(n)=af((alpha(nj)-mOC)/sOC4nf); 

pOD(n)=erf((alpha(nJ)-mOD)/sOD4nf); 

plA(n)=erf((mlA-alpha(nj))/slA4nf); 

plB(n)=erf((mlB-alpha(nj)yslB4nf); 

plC(n)=CTf((mlC-alpha(nJ))/slC,inf); 

plD(n)=etf((mlD-alpha(nj))/slD4nf); 

PBE(nj)=(pOA(n)+plA(n)+pOB(n)+plB(n)+pOC(n)+plC(n)+pOD(n)+plD(n))/16; 
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%  n»m+l; 


end  %Psloop 
% 

end  %Ckloop 

alsalphaO.l); 

a2«alpha(:^); 

a3salpha0.3); 

a43^pha(:,4); 

aSa=alphaO.S); 

a6salpha(:,6); 

a7salpha0,7); 

a8salpha(:,8); 

a9salpha0.9); 

alOBalpha(:,10); 

all«alpha(:.ll): 

pl*(PBE(:.l)); 

p2-(PBE(:^)); 

p3=(PBE(:3)); 

p4=(PBE(:.4)); 

p5=(PBE(:^)); 

p6=(PBE(:,6)); 

p7=(PBE(:.7)); 

p8=(PBE(:.8)); 

p9=(PBE(:.9)); 

pl0=PBE(:.10); 

pll=PBE(:.ll); 

PB=10*logl0(ps); 

plot(PB,al,'-'J>B^.'-'J>B,a3,’-’J>B.a4,’-’^B,a5;-’4>B,a6,’-'^B,a7, 

gridy«* 

%titie('case  1.3  (nu=.l,ddk«.3)'),.. 

xlabeK'Peak  Power  (dBW)'),ylabclCOptimal  Threshold’) 

s=[’Crosstalk=']; 

text(.15,.6,s,'sc') 

cl»0; 

$=['  ',num2str(cl)]; 
text(.2,.SS,s,'sc') 

[xs.ys]sdc2sc(PB(48).al(48)); 

polyline([.3,xs],[.56,ys],’-r',’sc’) 


c2«.35; 

ss[nuin2str(c2),'  dB']; 
text(.2,.S2,s,'sc') 
[xs,ys]»dc2sc(PB(48),a2(48)); 
polyline([.3,xs],[.53,ys],'-r','sc') 
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c3-32; 

s>s[nuin2str(c3),'  dB*]; 
textC2,.49,s,'sc') 
[xs,ys]»dc2sc(PB(48)^(48)); 
polyline(I.3.xs].[-5,ys],'-r’,’sc') 

c4»-29; 

sae[niiin2str(c4),'  dB*]; 
text(.2,.46,s,*sc') 
[xs,ys]sdc2sc(PB(48).a4(48)); 
polylinc([.3,xs],[.47.ys].’-r‘;sc’) 

cS»-26; 

ss[num2str(cS),'  dB']; 
text(.2,.43,s,'sc') 
[xs.ys]-dc2sc(PB(48),a5(48)); 
polyline([.3,xsM.44.ys].'-r’.’sc') 

c6»-23; 

s=[nuin2str(c6),'  dB']; 
text(.2,.4,s,'sc') 

[xs,ys]=dc2sc(PB(48).a6(48)); 

polylinc([.3^s],I.41.ys],'-r','sc') 


c7=-20; 

s=[nuin2str(c7),'  dB']; 
text(.2,.37,s,'sc') 
[xs.ys]=dc2sc(PB(48),a7(48)); 
polylinc([.3,xs],[.38,ys],'-r','sc') 

meta  t2chl3xa 


pause 


axis([-60-44-21  -5]) 

plot(PB,pl,'-'J>B.p2.'-'J>B.p3.'-'J>B.p4,'.'J>B.p5,'-'.PB,p6,'-'J>B.p7,'-')... 

semilogy.grid,.- 

%title('Perfonnancc  Curves:  case  1.3  (nu=.l,ddlc=.3)’),.. 

]dabel('Peak  Power  (dBW)'),ylabelCProbability  of  Bit  Error*) 


ss['Crosstalks']; 

text(.lS,.6^,'sc') 


cl»-20; 

ss[num2str(cl),'  dB']; 
text(.2,.55,s,'sc') 
[xs,ys]=dc2sc(PB(38).p7(38)); 
polyline([.3,xs],[.56,ys],'-r','sc') 


c2=-23; 

ss[nuin2str(c2),'  dB']; 
text(,2,.52,s,'sc') 
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(xs,ysl»dc2sc(PB(38),p6(38)); 

poIylinc(C.3^s],[.53,ys],’-r’,'sc') 

c3s-26; 

ss[nain2str(c3),'  dB*]; 
text(.2,.49.s.*sc') 
[xs.ys]=dc2sc(PB(38).p5(38)); 
polyline([.3,xs],[.5,ys],’-i:','sc') 


s«[num2str(c4),'  dB']; 
text(.2..46.s,'sc') 
(xs,ysl*sdc2sc(PB(38),p4(38)); 
polylinc([.3,xs],[.47,ys],'-r','sc’) 

c5*-32; 

s«[num2str(cS),'  dB']; 
tcxt(.2,.43.s,'sc') 
[xs,ys]*dc2sc(PB(38),p3(38)); 
polyline([,3,xs],[.44,ys],’-r','sc') 

c6a»-35; 

s«[nuni2str(c6),'  dB']; 
tcxt(.2,.4,s,'sc') 

[xs,ys]=dc2sc(PB(38),p2(38)); 

polylinc([.3,xs],[.41,ys],'-r','sc') 

c7=0; 

s=['  '4»um2str(c7)]; 
tcxt(.2,.37,s,'sc') 
(xs,ys]=dc2sc(PB(38),pl  (38)); 
polylinc((.3,xs],(.38,ys],'-r','sc') 
axis; 


meta  t2chl3xb 


APPENDIX  D  -  MATLAB  MODEL  OF  AN  OPTICAL  CHIP 
INTERCONNECT  FOR  ASYNCHRONOUS  TRANSMISSION 


This  model  is  fcv  the  case  of  2  adjacent  channels  with  (v,5j,52)  «  (0.1,0.3,0.3)  for 
asynchronous  transmission.  The  MATLAB  version  used  was  version  3.5k.  The  program 
itself  is  a  translation  of  equations  (14)  through  (25)  into  MATLAB  code. 


%  Asynchronous  Case:  Nos  lE-22 
%  ks2  implying  two  channels 

%  Case  1.3:  nus.l  and  delta  function  subscript  k  (delk)s.3 
%  th2asl3x.m 
%  20  July  93 

%  Known  constants 

Rs.5; 

T=1/(500E6); 

N0=lE-22; 

q=1.6E-19; 

Idk=10E-9; 

%  Coupling  constants  for  Ck^2 
ckl=0; 

ck2«3.16E-4;  %-35dB 
ck3»6.32E-4;  %-32dB 
ck4*1.26E-3;  %-29dB 
ck5-2.52E-3;  %-26dB 
ck5as5.011E-3:  %-23dB 
ckds.Ol;  %-20dB 
ck7=1.995E-2;  %-17  dB 
ck8s3.981E-2;  %-14dB 
ck9=7.943E-2;  %-ll  dB 
cklOs.1585;  %-8dB 

ck=[0  3.16E-4  6.32E-4  126E-3  2.52E-3  5.01  lE-3 .01  1.995E-2  3.981E-2  7.943E-2 
.1585); 

%  Patton  bit  values 

%bl=[.2182 .9759 .9759 .2182 .2182 .2182 .9759 .9759 .2182 .2182 .2182 .9759 
.9759 .9759 .2182 .9759); 

%b2=[.9759 .2182 .9759 .2182 .2182 .9759 .2182 .9759 .2182 .2182 .9759 .2182 
.9759 .2182 .9759 .9759); 

%fcfr  b32,b22,bl2,b02  for  patterns  A3,CJ>,EJ',G34  JJLL.M,N,04* 

%b3=[.9759 .9759 .9759 .9759 .2182 .2182 .2182 .2182 .2182 .9759 .2182 .2182 
.2182 .9759 .9759 .9759); 

%forb33,b23,bl3.b03 

%b4=[.9759 .9759 .9759 .9759 .9759 .9759 .9759 .9759 .2182 .2182 .2182 .2182 . 
2182 .2182 .2182 .2182); 
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%forb34,b24.bl4.b04 

blA»^182; 

blB».9759; 

blC«.9759; 

blD«.2182; 

blE».2182; 

blF«^182; 

blG-.9759; 

blH=.9759; 

bll=^182; 

blJ>.2182; 

blKs.2182; 

blL».9759; 

blM=.9759; 

blN=.9759; 

blCN.2182; 

blP=.9759; 

b2A=.9759; 

b2B».2182; 

b2C=.9759; 

b2D=.2182; 

b2E=.2182; 

b2F=.9759; 

b2G=.2182; 

b2H=.9759; 

b2I=.2182; 

b2J».2182; 

b2K=:.9759; 

b2M=.9759; 

b2N=.2182; 

b20=.9759; 

b2P=.9759; 

b3A=.9759; 

b3B=.9759; 

b3C=.9759; 

b3D=.9759; 

b3E=.2182; 

b3F=.2182; 

b3G=.2182; 

b3H=2182; 

b3I=.2182; 

b3J=.9759; 

b3K»^182; 

b3L=.2182; 

b3M=.2182; 

b3N=.9759; 

b30=.9759; 

b3P*.9759; 

b4A=.9759; 

b4B=.9759; 

b40.9759; 

b4D=.9759; 

b4E=.9759; 

b4F=.9759; 
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b4G-.9759; 

b4H».9759; 

b4I-^182; 

b4J«.2182; 

MK»^182; 

t'4L=.2182; 

o4M^.2182; 

MN=^182; 

b40.2182; 

b4P*.2182; 


b5A»blA'^2+b2A''2+b3A''2+b4A'^2; 

b5B=blB'^2+b2B^2+b3B^2+b4B''2; 

b5C=blC''2+b2C''2+b3C''2+b4C^2; 

b5D»blD^2+b2D^2+b3D^2+b4I>  >2; 

b5E*blE^2+b2E'^2+b3E'^2+b4E''2; 

b5F=blF^2+b2F''2+b3R2+b4F^2; 

b5G=blG''2+b2G''2+b3G''2+b4G''2; 

b5H=blH''2+b2H''2+b3H''2+b4H'^2; 

b5I=blH'^2+b2H''2+b3H''2+b4I''2; 

b5J»blJ''2+b2J''2+b3J''2+b4JA2; 

b5K=blKA2+b2K''2+b3K''2+b4K''2; 

b5M=blM''2+b2M''2+b3M'^2+b4M''2; 

b5N=blN''2+b2N'^2+b3N''2+b4N''2; 

b50=bl0^2+b20^2+b30^2+b40^2; 

b5P=blP^2+b2P''2+b3P^2+b4P'^2; 

bOO=..  .62\  %  approximate  zero 

bl0*.9759;  %  approximate  I 

nus.l; 

delks.3; 

ka=(pi*nu)/(4*pi^2*(nu''2+delk''2)); 

kb=l/(64*pi''6*(nu''2+<lelk''2)^3); 

kc=2*pi*nu*sin(2*pi*delk)+2*pi*delk*cos(2*pi*dclk); 

kd=(8*pi''2*nu*delk)*(2*pi*delk-(exp(-2*pi*nu))*kc); 

ke=(2*pi*delk)*sin(2*pi*delk)-(2*pi*nu)*cos(2*pi*delk); 

kf=4*pi''2*(delk''2-nu'^2)*(2*pi*nu-4*pi'^2*(nu^2+delk''2)+exp(-2*pi*nu)*ke); 

kska-kb*(kd+kf); 

n=l; 

j=i; 

m=l; 

%  •60dB<sPeak  Power<=-44dB 
%  Peak  power  is  A'^2/2  =  Ps 
for  ss=-60:.5:-44; 
ps(m>=l(yXss/10); 
m*m+l; 
end 

%  The  goal  here  is  to  solve  the  equation  for  the  optimal  threshold 
%  vs  p^  power 

alpha=zeTOsGength(ps),length(ck)); 
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PBEszeFOs(length(ps).length(ck)); 

for  js  1  :length(ck);  %  coupling  values  loop 
n=l; 

for  nal:length(ps);  %  peak  power  values 
Ps«ps(n); 

%  Constants 

qcnst=q/(Ps*T*R); 

vnoise=N0/(R''2*Ps'^2*T); 

vdark=(q*Idk)/(R'^2*Ps''2*T); 

%  Determine  the  sigma  for  each  of  the  0  bit  patterns 

s0A=sqrt(4*ck(j)*b0(y'2*k*b5A+qcnst*(b00^24ck(j)/2*b5A)+vnoise+vdark); 

s0B=sqn(4*ck(j)*b00^2*k*b5B+qcnst*(b00^2+ck(j)/2*b5B)+vnoise+vdark); 

sOC=sqrt(4*ck(j)*b00^2*k*b5C+qcnst*(bOO^?+ckO)/2*b5C)+vnoise+vdark); 

s0D=sqrt(4*ck0)*b00^2*k*b5D+qcnjt*{b00^2+ck0)^*b5D)+vnoise+vdark); 

s0Eassqrt(4*ck(j)*b00^2*k*b5E+qcnst*(b()(>^2+ck(j)^*b5E)+vnoise+vdark); 

s0F=sqn(4*ck^)*b00^2*k*b5F+qcnst*(b0(y^2+ck(j)/2*b5F)+vnoise+vdark); 

s0G=sqrt(4*ck(j)*b0Cy'2*k'''b5G+qcnst*(b(X)^2+ck(j)/2*b5G)+vnoise+vdark); 

s0H=sqrt(4*ck0)*b(Xy^2*k*b5H4qcnst*(b0(y'2+ck^)/2*b5H)+vnoise+vdark); 

s0I=sqrt(4*ck(j)*b0(y'2*k*b5I+qcnst*(b0(y'2+ck(j)/2*b5I)+vnoise+vdark); 

s0J=sqrt(4*ckQ'''b00^2*k*b5J+qcnst*(b00^2+ck0)^*b5J)+vnoise+vdark); 

s0K=sqn(4*ck(j)*b0(y'2*k*b5K+qcnst*(b00^2+ck(j)/2*b5K)+vnoise+vdark); 

sOL=sqrt(4*ck(j)*bO(y'2*k*b5L+qcnst»(bO(y'2+ck(j)^*b5L)+vnoise+vdaTk); 

s0M=sqrt(4*ck(j)*b00^2*k*b5M+qcnst*(b00^2+ck0)/2*b5M)+vnoise+vdark) 

s0N=sqrt(4*ck(j)*b00^2*k*b5N+qcnst*(b0(y^2+ck(j)^*b5N)+vnoise+vdark); 

s0O=sqrt(4*ck0)*b00^2*k*b5C)+qcnst*(b00'^2+ck(j)/2*b5O)+vnoise+vdark); 

s0P*sqn(4*ck(j)*b00^2*k*b5P+qcnst*(b00^2+ck(j)/2*b5P)+vnoise+vdark); 

%  Determine  the  variance  for  each  of  the  0  bit  patterns 


sOA2=sOA''2; 

sOB2=sOB''2; 

sOC2=sOC''2; 

sOD2=sOD'^2; 

sOE2=sOE''2; 

sOF2=sOF''2; 

sOG2=sOG^2; 

sOH2=sOH''2; 

s0I2=s0r2; 

sOJ2=sOJ''2; 

sOK2=sOK''2; 

sOL2=sOL''2; 

sOM2=sOM''2; 

sON2=sON''2; 

s002=s002; 

sOP2=sOP^2; 

%  Determine  the  sigma  for  each  of  the  1  bit  patterns 

slA=sqrt(4*ck(j)*bl(y'2*k*b5A+qcnst*(bl(y'2+ck(j)/2*b5A)+vnoise+vdark); 

slB=sqrt(4*ck^)*bl(y'2*k*b5B+qcnst*(bl(y^2+ck(j)/2*b5B)+vnoise+vdark); 
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slOsqrt(4*ckO)*bl(y'2*k*b5C4qcnst*(bl0^2+ck(j)/2*b5C)+vnoisc+vdark); 

slI>»sqrt(4*ck(j)*bl0^2*k*b5D4<icnsi*(bl(y'24ck(j)/2*b5D)+vnoise+v(iark); 

8lE»sqrt(4*ck(j)*bl(y'2*k*b5E+qcnst*(bl0^2+ck(iy2*b5E)+vnoisc+vdaik); 

slF»s<jrt(4*ck^)*bl(y'2*k*b5F+qcnst*^10^24ck(j)/2*b5F)+vnoise+vdaric); 

slG*sqrt(4*ck^)*bl0^2*k*b5G'Kicnst*(bl(y'2+ck(j)i/2*b5G)+vnoisc+vdark); 

slH=sqrt(4*ck(j)*bl0^2*k*b5H+qcnst*(bl(y'24ckOy2*b5H)+vnoisc+vdark); 

sll«sqrt(4*ck(j)*bl0^2*k*b5l4qcnst*(bl0^2+ck(j)fe*b51)+vnoisc4-vdaik); 

slJ»sqn(4*ck(j)*bl0^2*k*b5J4qcnst*(bl0^24ck(jy2*b5J)+vnoise+vdark); 

slK»sqrt(4*ckg)*bl(y'2*k*b5K+qcnst*(bl(y'2-K:ka)/2*b5K)+vnoise+vdark); 

slL»sqrt(4*ck(j)*bl0^2*k*b5L-Kicnst*(bl0^2+ck(jy2*b5LHvnoise+vdaik); 

slM«s^(4*ckg)*bl(y'2*k*b5M+qcnst*(bl0^24ck(jy2*b5M)+vnoisc+vdark); 

slN=s<]p(4*ckg)*bl(y'2*k*b5N+qcnst*(bl0^2+ckO)i^*b5N)+vnoisc+vdark); 

slO=sqit(4*ck(])*bl(y^2*k*b5O+qcnst*(bl(y^2+ck0)i^*b5O)+vnoise+vdark); 

slP=sqrt(4*ckg)*bl(y'2*k*b5P+qcnst*(bl0^2+ck(j)/2*b5P)+vnoise+vdark); 

%  Determine  the  variance  for  each  of  the  1  bit  patterns 


slA2=slA'^2; 

slB2=slB'^2; 

slC2»slC''2; 

slD2=slD^2; 

slE2=slE'^2; 

slF2=slF^2; 

slG2=slG''2; 

slH2=slH'^2; 

slI2=sir2; 

slJ2=slJ''2; 

slK2=slK'^2; 

slL2«slL^2; 

slM2=slM^2; 

slN2*slN'^2; 

8102=510^2; 

slP2=slP^2; 


%  Determine  the  mean  for  each  of  the  0  bit  patterns 

m0A=b00^2+ck(j)/2*b5A; 

m0B=b00^2+ckg)/2*b5B; 

mOC=bOO^2+ck0)/2*b5C; 

mOD=bOO^2+ck0)/2*b5D; 

mOE=bO(y'2+ck(j)/2*b5E; 

mOF=bO(y'2+ck(j)/2*b5F; 

mOG=b00^2+ck(j);2*b5G; 

mOH=bO(y'2+ckg)/2*b5H; 

m0I=b00^24ck(j)/2*b5I; 

m0J=b00^2+ck(j)/2*b5J; 

m0K=b00^2+ck(j)/2*b5K; 

mOL=bO(y'2+ckO)/2*b5L; 

m0M»b00^2+ckg)/2*b5M; 

m0N=b00^2+ck(j)/2*b5N; 

m00=b00^2+ck(j)/2*b50; 

m0P=b00^2+ck(j)/2*b5P; 
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%  Determine  the  mean  for  each  of  the  1  bit  patterns 


mlA»bl0^2+ck(j)/2*b5A; 

mlB=bl0^24ckO)/2*b5B; 

mlC-bl0^24ckO)/2*b5Q 

mlD»bl0^2+ckO)/2*b5D; 

mlE=bl(y'2+ckO)/2*b5E; 

mlF-bl0^2+ck(j)/2*b5F; 

mlG«bl0^2+ck(j)/2*b5G; 

mlH«bl0^2+ck(j)/2*b5H; 

mll»bl0^24ck(j)/2*b51; 

mlJ»bl0^2+ck(j)/2*b5J; 

mlK=bl0^2-K;k(j)/2*b5K; 

mlL=bl0^2+ck(j)/2*b5L; 

mlM=bl0^24ck(j)/2*b5M; 

mlN=bl0^2+ck(j)/2*b5N; 

mlO=bl0^2+ck(j)/2*b5O; 

mlP=bl0^2-K:k(j)/2*b5P; 

%Solve  the  equation  F(a)=0  to  optimize  a,  the  threshold.  F(a)=0  comes 

%  from  (dP/da)sO 

a— 0*  0009*1' 

FA=(l/s0A)*exp(-((a-m0A).''2)./(2*s0A2))-(l/slA)*exp(-((mlAa).''2)/(2*slA2)); 

FB=(l/s0B)*cxp(-((a-m0B).A2)./(2*s0B2))-(l/slB)*exp(-((mlB-a).''2)/(2*slB2)); 

FC=(l/s0C)*cxp(-((a-m0C).''2)./(2*s0C2))-(l/slC)*cxp(-((mlC-a).^2)/(2*slC2)); 

FD=(l/s0D)*exp(-((a-m0D).''2)./(2*s0D2))-(l/slD)*cxp(-((mlD-a).''2)/(2*slD2)); 

FE=(l/s0E)*exp(-((a-m0E).^2)./(2*s0E2))-(l/slE)*exp(-((mlE-a).''2)/(2*slE2)); 

FF=(l/sOF)*exp(-((a-mOF).''2)y(2*sOF2))-(l/slF)*cxp(.((mlF-a).^2)/(2*slF2)); 

FG=(l/s0G)*cxp(-((a-m0G).^2)./(2*s0G2))-(l/slG)*cxp(-((mlG-a).'^2)/(2*slG2)); 

FH=(l/s0H)*cxp(-((a-m0H).''2)./(2*s0H2))-(l/slH)*cxp(-((mlH-a).''2)/(2*slH2)); 

H=(l/s0I)*cxp(-((a-m0I).''2)./(2*s0I2))-(l/slI)*cxp(-((mlI.a).'^2)/(2*slI2)); 

FJ=(l/s0J)*cxp(-((a-m0J).^2)./(2*s0J2)Hl/slJ)*cxp(.((mlJ-a).'^2)/(2*slJ2)); 

FK=(l/s0K)*cxp(-((a.m0K).''2)./(2*s0K2))-(l/slK)*cxp(-((mlK-a).''2)/(2*slK2)); 

FL=(l/s0L)*cxp(-((a-m0L).''2)./(2*s0L2))-(l/slL)*cxp(-((mlL-a).^2)/(2*slL2)); 

FM=(l/s0M)*exp(-((a-m0M).''2)./(2*s0M2))-(l/slM)*cxp(-((mlM- 

a).^2)/(2*slM2)); 

FN=(l/s0N)*cxp(-((a-m0N).'^2)./(2*s0N2))-(l/slN)*cxp(-((mlN-a).’^2)/(2*slN2));’ 

F0=(l/s00)*cxp(-((a-m00).'^2)./(2*s002))-(l/sl0)*cxp(-((ml0-a).^2)/(2*sl02)); 

FP=(l/s0P)*cxp(-((a-m0P).'^2)y(2*s0P2))-(l/slP)*cxp(-((mlP-a).'^2)/(2*slP2)); 

F=FA+FB+FC+FD+FE+IT+FG+FH+FI+FJ+FK+FL+FM+FN+FO+FPj 

[Y,I]=min(abs(F)); 


alpha(nj)=a(l); 

aa=alpha(nj); 

%  Determine  the  probability  of  bit  error  for  each  optimal  threshold 
%  over  all  bit  patterns  A,B.CJD 
pOA(n)=crf((alpha(no)-mOA)/sOA,inO; 
p0B(n)=erf((alpha(nj)-m0B)/s0B,inO; 
pOC(n)=erf((alpha(nJ)-mOC)/sOC,inf); 
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pOD(n)*erf((alpha(nJ)-mOD)/sOD.inf); 

pOE(n)aerf((alpha(nJ)-mO£)/sOEjnf); 

pOF(n)Beif((alpha(nj>mOF)/sOF,inf): 

pOG(n)a«rf((alpha(nj)-mOG)/sOG.if^; 

pOH(n)«erf((alpha(nJ)'n^H)/sOH.inf); 

pOI(n)aeif((dpna(nJ)-mOI)/^I4nf); 

pQJ(n)aerf((alpha(no>inOJ)/sOJ4nf): 

pOK(n)-erf((a]pha(no)’mOKVsOK4nf): 

pOL(n)««rf((alpha(nJ)-mOL)/sOLinf); 

pOM(n)serf((alpha(nj)-ii]OM)/sOM.iiif); 

pON(n)»crf((alpha(nj)-mON)/sON,inf); 

p0C>(n)s«rf((alpha(nJ)-m0OVs0O4nf): 

pOP(n)seif((alpha(nj>*mOP)i/sOP,inf); 

plA(n)aerf((inlA-alpha(nJ))/slA.inf); 

plB(n)=crf((mlB-alpha(nJ)yslB,inf); 

plC(n)s=crf((mlC-aIpha(nj))/slC,inf); 

plD(n)serf((mlD-alpha(nJ))/slD,inO; 

plE(n)serf((n)lE*alpha(nJ))/slE,inO; 

plF(n)=crf((mlF-alpha(nj))^lF,inf); 

plG(n)s«rf((mlG*alpha(nJ))/slG,ii^; 

plH(n)=erf((inlH-alpha(nj)VslH.inf); 

pll(n)serf((mll-alpha(nj))/sll4n0; 
plJ(n)s=erf((mlJ-alpha(nj))/slJ,inO; 
p  1  K(n)=erf((ml  K-dpha(n  J))/s  1  K,inf); 
plL(n)=crf((mlL-alpha(nj))/slL,inO; 
plM(n)serf((inl  M-^pha(n  J))/s  1  M,inO; 
plN(n)=crf((mlN-alpha(no))/slN,inf); 
p  1 0(n)=af((ml  0-alpha(n  J)ysl  0,inf); 
plP(n)=crf((mlP-alDha(nj))/slP,inO; 

pA(n)=pOA(n)+plA(n); 

pB(n)=pOB(n)+plB(n); 

pC(n)=pOC(n)+plC(n); 

pD(n)=pOD(n)+p  lD(n); 

pE(n)=^E(n)+plE(n); 

pF(n)=pOF(n)+plF(n); 

pG(n)=pOG(n)+plG(n); 

pH(n)=pOH(n)+plH(n); 

pl(n)=i^)l(n)+pll(n); 

pJ(n>:pOJ(n)+plJ(n); 

pK(n)=pOK(n)+plK(n); 

pL(n)=^L(n)+plL(n); 

pM(n)=pOM(n)+p  1  M(n); 

pN(n)=^N(n)+plN(n); 

pO(n)=pOO(n)+p  10(n); 

pP(n)=pOP(n)+plP(n); 


PBE(nj)=(pA(n)+pB(n)+pC(n)+pD(n)+pE(n)+pF(n)+pG(n)+pH(n)+pI(n)+pJ(n)+pK 

(n)+pL(n)+pM(n)+pN(n)+pO(n)+pP(n))/64; 


end  %Psloop 
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end  %C!k  loop 


al«alpha(:,l); 

a2>«lpha(:^); 

a3«alpha(:,3); 

a4a<alpha(:,4); 

aS»alpha(:4); 

a6>alpha(:,6); 

a7a>alpha(:,7); 

a8=alpha(^8); 

a9salphaC.9); 

al0salpha(:,10): 

all«alpha(;,ll); 

pmPBE(:.l)); 

p2=(PBE(.%2)); 

p3=(PBE(;.3)); 

p4=(PBE(:.4)); 

p5=(PBE(:,5)); 

p6=(PBE(:.6)); 

p7=(PBE(:.7)); 

p8=(PBE(:.8)); 

p9=(PBE(:.9)); 

pl0=PBE(:.10); 

pll=PBE(:.ll); 


PB=10*logl0(ps); 

plot(PB.ai;-’,PB,a2;-’J»B.a3;-'.PB.a4yj>B.a5;-’J>B,a6;.'JPB,a7;.’),.. 

giidf** 

%title('casc  1.3,nu*.l,delk=.3:  Fovir  channels-Asynchronous 
xlabeU’Pcak  Powcr(dBW)'),ylabel('Optimal  Threshold’) 


s=['Crosstalk=']; 

tcxt(.15,.6,s,'sc') 

cl=0; 

ss['  '4ium2str(cl)]; 
tcxt(.2,.55,s,'sc’) 
[xs.ys]=dc2sc(PB(19)^l(19)); 
polyline([.3,xs],[.56,ys],'-r','sc') 

c2=-35; 

s=[num2str(c2),’  dB']; 
text(.2,.S2,s,'sc') 
[xs,ys]=dc2sc(PB(19)^(19)); 
polyline(I.3,xs],[-53,ys],'-r’  ,'sc') 


c3=-32; 

s=[nuin2str(c3),'  dB’]; 
tcxt(.2,.49,s,'sc') 
[xs.ys]=dc2sc(PB(19),a3(19)); 
polyline([.3,xs],1.50,ys],'-r','sc’) 
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c4^29; 

s»[niini2str(c4),'  dB']; 
text(.2..46.s,'sc') 

[xs,y  s]a>dc2sc(PB(  19),a4(  1 9)); 
polylinc([.3,xs],[.47,ys],’-r',*sc') 

C5-26; 

s>:[num2str(cSV  dB']; 
text(.2,.43,s.'sc') 
[xs,ys]aKlc2sc(PB(19),a5(19)); 
polyliiic([.3,xs],[.44,ys],’-r',*sc') 

c6=-23; 

ss[nuni2str(c6),'  dB']; 
text(.2,.40,s.'sc') 
[xs,ys]»dc2sc(PB(19)^6(19)); 
polyliiie([.3,xs].[.4 1  ,ys],'-r','sc') 


c7=-20; 

ss[nuin2str(c7),'  dB']; 
tcxt(.2,.37,s,'sc') 

[xs,ys]=dc2sc(PB(19).a7(19)); 

polyline([.3,xs],[.38,ys],'-r','sc') 

metat2asl3xa 

pause 

axis([-60-44-21  -5]) 

plot(PB.pl.'-'JPB.p2.'-'4>B.p3;-'JPB.p4,'-'4>B.p5.'-'4>B,p6,'^'.PB,p7.'.'),.. 

semilogy.grid,.. 

%title('case  1.3,nu=.l,delk=.3:  Four  channcls-Asynchronous'),.. 
xlabeK'Peak  Power(dBW)'),ylabcl('Probability  of  Bit  Error') 


s=['Crosstalk=']; 

text(.15,.6,s,'sc') 


cl=-20; 

s=[nuixi2str(cl),'  dB']; 
tcxt(.2,.55,s,'sc') 
[xs,ys]=dc2sc(PB(15).p7(15)); 
polylinc([.3,xs],[.56,ys],'-r’,'sc') 

c2=-23; 

ss[nuin2str(c2),'  dB']; 
text(.2,.52,s,'sc') 
[xs,ys]»dc2sc(PB(15),p6(15)); 
polylinc([.3^],[.53,ys],'-r','sc') 

c3=-26; 

ss[num2str(c3),'  dB']; 
tcxt(.2,.49,s,'sc') 
[xs,ys]=dc2sc(PB(15).p5(15)); 
polyline([.3,xs],[.50,ys],'-r','sc') 
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c4»-29; 

s>B[nuiii2str(c4),*  dB']; 
text(^,.46,s,'sc') 

[xs  ,ys]«dc2sc(PB(  1 5),p4(  1 5)); 
polylinc([.3pts],[-47,ys],'-r'.’sc') 

C5-32; 

s«[num2str(c5),'  dB']; 
text(.2..43,s,'sc') 
[xs.ys]=dc2sc(PB(15).p3(15)); 
polyline([.3^s],[.44,ys],'-r’,'sc') 

C6-35; 

ss[nuin2str(c6),'  dB']; 
text(.2,.40,s,'sc') 
[xs,ys]=dc2sc(PB(15).p2(15)); 
polylinc([.3,xs],[.4 1  .ys],'-r’,'sc') 

c7=0; 

ss['  '^uGQ2str(c7)]; 

tcxt(.2,.37,s,'sc') 

[xs,ys]=dc2sc(PB(15),pl(15)); 

polyline([.3,xs],[.38,ys],'-r’,'sc') 

axis; 

tneta  t2asl3xb 
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APPENDIX  E  -  MATLAB  MODEL  FOR  THREE  ADJACENT 
CHANNELS  WITH  SYNCHRONOUS  TRANSMISSION 


This  model  is  for  the  case  of  three  adjacent  channels  with  (v,5i.52.53)  « 
(0.1,0.3,0.3,0.3)  for  synchronous  transmission.  The  MATLAB  version  used  was  version 
3.Sk.  The  program  itself  is  a  translation  of  equations  (1)  through  (13)  into  MATLAB  code 
for  k«3. 


%  Synchronous  Case 
%  ks3  implying  three  channels 

%  Case  1.3:  nu^.l  and  delta  function  subscript  k  (delk)s.3 
%  th3chl3.m 
%  1  Aug  93 

%  Thermal  Noisc=lE-22  A'^Z/Hz 

%  Known  constants 

R=.5; 

T=1/(500E6); 

N0=lE-22; 

q=1.6E-19; 

Idk=10E-9; 

%  Coupling  constants  for  Ck^2 
ckl=0; 

ck2«3.16E-4;  %-35dB 
ck3=6.32E-4;  %-32dB 
ck4=1.26E-3;  %-29dB 
ck5=2.52E-3;  %-26dB 
ck5a=5.011E-3;  %-23  dB 
ck6=.01;  %-20dB 

ck7=1.995E-2;  %-17  dB 
ck8=3.981E-2;  %-14dB 
ck9=7.943E-2;  %-ll  dB 
ckl0=.1585;  %-8dB 

ck=[0  3.16E-4  6.32E-4  1.26E-3  2.52E-3  5.012E-3  .01  1.995E-2  3.981E-2 
7.943E-2  .1585]; 

%  Pattern  bit  values 

%bl«[.2182  .9759  .9759  .2182  .2182  .2182  .9759  .9759]; 

%for  b21,bll,b01  for  patterns  A,B.CJD3J^.GJ1 
%b2=[.9759  .2182  .9759  .2182  .2182  .9759  .2182  .9759]; 

%for  b22,bl2,b02  for  patterns  A,B,C4)3J,G3 
%b3=[.9759  .9759  .9759  .9759  .2182  .2182  .2182  .2182]; 

%for  b23,bl3,b03 
blA».2182; 
blB=.9759; 
blC».9759; 
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bll>»Jil82; 

blE«.2182; 

blF*.2182; 

blG«.9759; 

blHs.9759; 

b2A*.9759; 

b2B».2182; 

b2C».9759; 

b2D».2182; 

b2E«.2182; 

b2F».9759; 

b2G«.2182; 

b2H*.9759; 

b3A=.9759; 

b3Bs.9759; 

b3C*.9759; 

b3Ds.9759; 

b3E».2182; 

b3F».2182; 

b3G«.2182; 

b3H=.2182; 

b4A=blA''2+b2A''2+b3A''2; 

b4B*blB'^2+b2B'^2+b3B'^2; 

b4C*blC'^2+b2C'^2+b3C''2; 

b4D=bliy'2+b2D^2+b3D^2; 

b4E*blE'^2+b2E'^2+b3E^2; 

b4F=blFA2+b2F''2+b3F'^2; 

b4G*blG''2+b2G'^2+b3G''2; 

b4H«bIH'‘2+b2H^2+b3H''2; 

bOO=.2182;  %  approximate  zero 
blO=.9759;  %  ai^roximate  1 

%  Detennine  k  constant  for  the  variance  of  X 

nus.l; 

delk«.3; 

ka=(2*pi*nu)/(4*pi''2*(nu''2+delk''2)); 

kb=l/(16*piM*(nu^2+delk''2r2); 

kc=(2*pi*nu*exp(-2*pi*nu))*(2*pi*delk*sin(2*pi*delk)- 

2*pi*nu*cos(2*pi'''deUc)); 

kd=<2*pi*delk*exp(-2*pi*nu))*(2*pi*(delk*cos(2*pi*delk)+nu*sin(2*pi*dclk))); 

ke=4*pi^2*(nu''2-delk^2); 

kska-kb^Ckc+kd+ke); 

n=l; 

j=i; 

m=:l; 


%  •60dB<BPeak  Power<='44dB 
%  Peak  power  is  A''2/2  =  Ps 
for  ss=-60:.2:-44; 
ps(m)=l(>Xss/10); 
iiis=m+l; 
end 

%  The  goal  here  is  to  solve  the  equation  for  the  optimal  threshold 
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%  vs  peak  power 


atoha"zeros(length(ps)4ength(ck)); 

PBEszeiosOength(ps)Jen^(ck)); 

for  j«l  :length(ck);  %  coupling  values  loop 
n»l; 

for  n»l  dength(ps);  %  peak  ppower  values 
I*s*ps(n); 

%  Constants 

qcnst^(Ps*T*R.); 

vnoise»N0/(R'^2*Ps''2*D; 

vdaik»(q*Idk)/(R'^2*Ps'^2*T); 

%  Deteitnine  the  sigma  for  each  ctf  the  0  bit  patterns 

sOA*sqrt(4*ck(j)*bO(y'2*k*b4A+qcnst*(bO(y'2-K;k(j)*b4A)+vnoise+vdaik); 

s0B=sqTt(4*ck(j)*b00^2*k*b4B+qcnst*(b00'^2+ck(j)*b4B)+vnoisc+vdark); 

s0C*sqrt(4*ck0)*b00^2*k*b4C4qcnst*(b00^2+ck(i)*b4C)+vnoise+vdark); 

sOD=s<pt(4*ckO)*bOO'^2*k*b4I>+qcnst*{bOO''2+ckO)*b4D)+vnoise+vdark); 

sOE»sqrt(4*ck(j)*bO(y'2*k*b4E+qcnst*(bO(y^2+ck(j)*b4E)+vnoise+vdark); 

sOF=s(p(4*ckO)*bO(y^2*k*b4F+qcnst*(bOO'^24ckO)*b4F)+vnoise+vdark); 

s0G*sqrt(4*ck0)*b00^2*k*b4G-Kicnst*(b0Cy^2+ck(j)*b4G)+vnoise+vdark); 

s0H»sqrt(4*ck(j)*b(Xy^2*k*54H+qcnst*(b00'^24ck0)*b4H)+vnoise+vdark); 

%  Detemiine  the  variance  for  each  of  the  0  bit  patterns 

sOA2*sOA''2; 

sOB2=sOB'^2; 

sOC2«sOC'^2; 

sOD2»sOD^2; 

sOE2*sOE''2; 

sOF2=sOF''2; 

sOG2»sOG'^2; 

sOH2=sOH^2; 


%  Detemiine  the  sigma  for  each  of  the  1  bit  patterns 

slA=sqrt(4*ck(j)*bl(y'2*k*b4A+qcnst*(bl0^2+ck(j)*b4A)+vnoise+vdaA); 

slB«sqrt(4*ckO)*bl(K'2*k*b4B4qcnst*(blO''2+ck(j)*b4B)+vnoise+vdaik); 

slC«sqrt(4*ckO)*bl0^2*k*b4C-Kicnst*(bl0^2+ck(j)*b4C)+vnoise+vdaik); 

slD*sqrt(4*ck(j)*bl0^2*k*b4D+qcnst*(bl0^2+ck(j)*b4D)+vnoise+vdaik); 

slE=s^(4*ck(j)*bl(y'2*k*b4E4qcnst*(bl0^2+ck(j)*b4E)+vnoise+vdark); 

slF»sqrt(4*ck(j)*blCX'2*k*b4F+qcnst*(bl0^2-K;k(j)*b4F)+vnoise+vdark); 

slG«s^(4*ckQ)*bl0'^2*k*b4G4qcnst*(bl(y'2+ck(j)*b4G)+vnoise+vdark); 

slH=sqrt(4*ckO)*bl(y'2*k*b4H+qcnst*(bl0^2+ckO)*b4H)+vnoise+v(lark); 

%  Determine  the  variance  for  each  of  the  1  bit  patterns 


slA2=slA'^2 

slB2:«slB^2 

slC2»slC''2 
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slD2»slD^2; 

slE2-slE^2; 

slF2»slF''2; 

slG2»slG^2: 

slH2«slH^2; 


%  Detennine  die  mean  for  each  of  the  0  bit  patterns 
m0A»b00^2+ck(j)*b4A; 
m0B=b00^2+ck(j)*b4B; 
m0Ob00^2+ck(j)*b4C 
m0D-b00^2-H:k(j)*b4D; 
m0E=b00^2+ck(j)*b4E; 
mOF»bO(y'24ckO)*b4F; 
m0G»b00^2+ck(j)*b4G; 
m0H=b00A2+ck(j)*b4H; 


%  Detennine  the  mean  for  each  of  the  1  bit  patterns 

mlA=bl(y^2+ck(j)*b4A; 

mlB=bl0^2+ck(j)*b4B; 

mlC=bl0^2+ck(j)*b4C; 

mlD=bl0^2+ckO)*b4D; 

mlE=bl0^2+ck(j)*b4E; 

mlF=bl(y'2-K:k(j)*b4F; 

mlG=bl(y'2+ck(j)*b4G; 

mIH=:bl0^24ck(j)*b4H; 

%Solve  the  equation  F(a)=0  to  optimize  a,  the  threshold.  F(a)=0  comes 
%  from  (dP/da)=0 

a=0:.001:l; 

FA=(l/s0A)*exp(-((a-m0A).''2)./(2*s0A2))-(l/slA)*exp(-((mlA-a).^2)/(2*slA2)); 

FB=(l/s0B)*exp(-((a-m0B).^2)./(2*s0B2))-(l/slB)*exp(-((mlB-a).''2)/(2*slB2)); 

FC=(l/s0C)*exp(-((a-m0C).''2)./(2*s0C2))-(l/slC)*exp(-((mlC-a).^2)/(2*slC2)); 

FD=(l/s0D)*exp(-((a-m0D).''2)./(2*s0D2)Hl/slD)*exp(-((mlD-a).'^2)/(2*slD2)); 

FE=(l/s0E)*exp(-((a-m0E).^2)./(2*s0E2))-(l/slE)*exp(-((mlE-a).''2y(2*slE2)); 

FF=(l/s0F)*exp(-((a-m0F).''2)./(2*s0F2))-(l/slF)*exp(-((mlF.a).^2)/(2*slF2)); 

FG»(l/s0G)*exp(-((a-m0G).''2)./(2*s0G2))-(l/slG)*exp(-((mlG-a).''2)/(2*slG2)); 

FH=(l/s0H)*exp(-((a-m0H).''2)./(2*s0H2))-(l/slH)*exp(-((mlH-a).''2)/(2*slH2)); 

F*FA+FB+FC+ro+FE+FF+FG+FH; 

[Y,I]=min(abs(F)); 

alpha(nj>4(l); 

aasalpha(nj); 

%  Detennine  Ae  probability  of  bit  enor  for  each  optimal  threshold 
%  over  all  bit  patterns  A,B>C 
pOA(n)=crf((alpha(nj)-mOA)/sOA4nO; 
pOB(n)*erf((alpha(nj)-mOB)/sOB4nf); 
pOC(n)=erf((alpha(no)-mOC)/sOC,inf); 
pOD(n)=erf((alpha(nj)-mOD)/sOD,inf); 
pOE(n)®eif((alpha(nj)-mOE)/sOE,inf); 


71 


pOF(ii)«erf((alpha(nJ)-mOFVsOF4nf); 

pOG(n)«eif((alpha(nJ)-mOG)/sOG4i^; 

pOH(n)«eif((alpha(iiJ)-mOH)/sOH4nf): 

plA(n)-eif((iDlA-alpha(nJ)yslAanf); 

plB(n)»eif((inlB-alpha(nJ)VslB4nf); 

plC(n)«erf((mlC-alpha(nJ)VslC4iif): 

plD(n>*erf((mlD-alplui(nJ)VslD4nf): 

plE(n)aeif((mlE-alpha(nJ)VslE4nf); 

plF(n)aerf((mlF-alpha(nj)VslF4nf); 

plG(n)seif((mlG*alpha(nJ))/slG4iii): 

plH(n)serf((mlH-alpha(nJ))/slH4nf); 


PBE(nj)»(pOA(n)+plA(n)+pOB(n)+plB(n)+pOC(n)+plC(n)+pOD(n)+plD(n) 

: .  : : .  :  .“'':n)+piG(n)+pOH(n)+piH(n))/32; 


end  %Psloop 

end  %Ck  loop 

alsalpha(:,l); 

a2salphaO»2V, 

a3=^phaC.3); 

a4=alpha(:,4); 

aS3^pha0.5); 

a6=alpha(:,6); 

a7=alpha(:,7); 

a8salpha(:,8); 

a9=alpha(:,9); 

al0salpha(:,10); 

all=alpha(:,ll): 

pl=(PBE(:.l)); 

p2=(PBE(:.2)); 

p3=(PBE(:.3)); 

p4=(PBE(:.4)); 

p5=(PBE(:.5)); 

p6*(PBE(:.6)); 

p7=(PBE(:.7)); 

p8=(PBE(:.8)); 

p9=(PBE(:.9)); 

plO=PBE(:,10); 

pll«PBE(:,ll); 

PB=10*logl0(ps); 

plot(PB.ai;-',PB.a2;-’.PB,a3,’-\PB,a4;-’,PB.a5,’-’,PB,a6.’-’,PB,a7 

grid... 

%titleCcase  1.3,nu=.l,delk=.3:  Three  channels’),.. 
xlabelC'Peak  Power(dBW)'),ylabel(’C)ptinial  Threshold’) 


s=[’Crosstalk=’]; 

text(.15,.6,s,’sc’) 
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cl-0; 

s^r  '4iuin2str(cl)]; 
text(.2..5S,s,'sc') 
[xs,ys]adc2sc(PB(48),a  1  (48)); 
polylinc([.3,xs],[.56,ys],'-r’,'sc') 

c2^35; 

s«[num2str(c2),'  dB']; 
text(.2..S2,s,'sc*) 
Ixs.ys]*=dc2sc(PB(48),a2(48)); 
polylinc([.3,xs],[.53,ys],’-r’,'sc') 

c3-32; 

s»[nua^str(c3),'  dB']; 
tcxt(.2,.49,s,'sc’) 
[xs,ys]sdc2sc(PB(48),a3(48)): 
polyline((.3,xs],(.5,ys].'-r','sc') 

c4=-'^; 

s>:[num2str(c4),'  dB']; 
text(.2..46,s.'sc') 
[xs,ys]=dc2sc(PB(48),a4(48)); 
polylinc((.3,xs],[.47,ys],'-r','sc') 


c5»-26; 

s=[num2str(c5),'  dB']; 
tcxt(.2,.43,s,'sc') 

Ixs.ys]=dc2sc(PB(48),a5(48)); 

polyline((.3,xs],[.44,ys],'-r','sc') 

c6=-23; 

s»[nuin2str(c6),'  dB']; 
text(.2,.4,s,'sc') 

[xs,ys]=dc2sc(PB(48),a6(48)); 

polyline([.3,xs],[.41,ys],'-r','sc') 

c7=-20; 

s=[num2str(c7),'  dB']; 
tcxt(.2,.37,s,'sc') 

[xs,ys]sdc2sc(PB(48),a7(48)); 
polylinc([.3,xs],[.38,ys],'-r','sc') 
meta  t3chl3xa 

pause 

axis([-60-44-21  -5]) 

plot(PB,pl,'-'.PB.p2,'-’,PB.p3,’-'.PB,p4;-'.PB.p5.'-’,PB,p6,'-'.PB.p7,'-'),.. 

seinilogy,grid,.. 

%title('case  1.0,nu=.l,delk=0:  Three  channels'),.. 
xlabelCPeak  Power(dBW)'),ylabelCProbability  of  Bit  Error') 


s=['Crosstalk=']; 

text(.15,.6,s,'sc') 
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cl-20; 

s»[num2str(cl),'  dB']; 
text(.2,.55,s,'sc’) 
[xs,ys]=dc2sc(PB(38),p7(38)); 
polylinc([.3,xsl,[.56,ys],'-r’,'sc') 

c2*-23; 

ss[nuiii2str(c2),'  dB']; 
tcxt(.2,.52,s,’sc') 
[xs,ys]*=dc2sc(PB(38),p6(38)); 
polylinc([.3,xs],[.53,ys],'-r','sc‘) 

c3»-26; 

ss[nuD^str(c3),'  dB']; 
tcxt(.2,.49,s,'sc') 
[xs,ys]=dc2sc(PB(38),p5(38)); 
polyline([.3,xs],[.5,ys],'-r','sc') 

s=[nuin2str(c4),'  dB']; 
text(.i,.46,s,’sc') 
[xs,ys]=dc2sc(PB(38).p4(38)); 
polylinc([.3,xsl,[.47,ys],'-r','sc’) 


c5=-32; 

s=[nuin2str(c5),'  dB']; 
tcxt(.2,.43,s,'sc') 
(xs,ys]=dc2sc(PB(38),p3(38)); 
poIyIine({.3,xs],[.44,ys],'~r','sc') 

c6=-35; 

s=[nuin2str(c6),'  dB']; 
text(.2,.40,s,'sc') 
[xs,ys]=dc2sc(PB(38),p2(38)); 
polylinc([.3,xs],[.41,ysl,'-r','sc') 

c7=0; 

s=['  ',nuin2str(c7)]; 
text(.2,.37,s,'sc') 
[xs,ys]=dc2sc(PB(38),pl(38)); 
polyline([.3,xs],[.38,ys],'-r','sc') 

axis; 

meta  t3chl3xb 
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APPENDIX  F  -  MATLAB  MODEL  FOR  FOUR  ADJACENT 
CHANNELS  FOR  SYNCHRONOUS  TRANSMISSION 


This  mcxiel  is  for  the  case  of  four  adjacent  channels  with  (v,5i, 52,53.54)  » 
(0.1,0.3,0.3,0.3,0.3)  for  synchronous  transmission.  The  MATLAB  version  used  was 
version  3.Sk.  The  program  itself  is  a  translation  of  equations  (1)  through  (13)  into 
MATLAB  code  for  k=4. 

%  Synchronous  Case 
%  ks4  implying  four  channels 

%  Case  1.3:  nu=.l  and  delta  function  subscript  k  (delk)=.3 
%  th4chl3x.m 
%  22  Aug  93 

%  Theim^  Noise:  lE-22  A'^2/Hz 

%  Known  constants 

R=.5; 

T=1/(500E6); 

N0=lE-22; 

q=1.6E-19; 

Idk=10E-9; 

%  Coupling  constants  for  Ck'^2 
ckl*0; 

ck2=3.16E-4;  %-35dB 
ck3=6.32E-4;  %-32dB 
ck4=1.26E-3;  %-29dB 
ck5=2.52E-3;  %-26dB 
ck5a=5.011E-3;  %-23dB 
ck6=.01;  %-20dB 

ck7=1.995E-2;  %-17  dB 
ck8=3.981E-2;  %-14  dB 
ck9=7.943E-2:  %-ll  dB 
cklO=.1585;  %-8dB 

ck=[0  3.16E-4  6.32E-4  1.26E-3  2.52E-3  5.01 1E.3  .01  1.995E-2  3.981E-2 
7.943E-2  .1585]; 

%  Pattern  bit  values 

%bl=1.2182  .9759  .9759  .2182  .2182  .2182  .9759  .9759  .2182  .2182  .2182 
.9759  .9759  .9759  .2182  .9759]; 

%for  b31,b21,bll,b01  for  patterns  A,B.CJ>,E,F,GJi,U,K,L,M,N.OJ» 
%b2=[.9759  .2182  .9759  .2182  .2182  .9759  .2182  .9759  .2182  .2182  .9759 
.2182  .9759  .2182  .9759  .9759]; 

%for  b32,b22.bl2,b02  for  patterns  A,B,CjDJE,F,GJl,I,J,KL,M,N,0,P 
%b3=[.9759  .9759  .9759  .9759  .2182  .2182  .2182  .2182  .2182  .9759  .2182 
.2182  .2182  .9759  .9759  .9759]; 

%for  b33,b23,bl3,b03 
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%b4»[.9759  .9759  .9759  .9759  .9759  .9759  .9759  .9759  .2182  .2182  .2182 
.2182  .2182  .2182  .2182  .2182]; 

%for  b34,b24,bl4,b04 

blAs.2182; 

blB=.9759; 

blC=.9759; 

bID».2182; 

blE=.2182; 

blF=.2182; 

blG».9759; 

blH=.9759; 

bll=.2182; 

blJ=.2182; 

blK=.2182; 

blL=.9759; 

blM=.9759; 

blN=.9759; 

blO.2182; 

blP=.9759; 

b2A=.9759; 

b2B=.2182; 

b2C=.9759; 

b2D=.2182; 

b2E=.2182; 

b2F=.9759; 

b2G=.2182; 

b2H=.9759; 

b2I=.2182; 

b2J=.2182; 

b2K=.9759; 

b2L=.2182; 

b2M=.9759; 

b2N=.2182; 

b20=.9759; 

b2P=.9759; 

b3A=.9759; 

b3B=.9759; 

b3C=.9759; 

b3D=.9759; 

b3E=.2182; 

b3F=.2182; 

b3G=.2182; 

b3H=.2182; 

b3I=.2182; 

b3J=.9759; 

b3K=.2182; 

b3L=.2182; 

b3M=.2182; 

b3N=.9759; 

b30.9759; 

b3P=.9759; 

b4A=.9759; 

b4B=.9759; 

b4C=.9759; 
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b4D=.9759; 

b4E«.9759; 

b4F».9759; 

b4G«.9759; 

b4H»,9759; 

b41».2182; 

b4J»^182; 

b4K«.2182: 

b4L«.2182; 

b4M»^182; 

b4N=.2182; 

b40.2182; 

b4P=.2182; 


b5A=blA'^2+b2A^2+b3A'^2+b4A'^2; 

b5B=blB'^2+b2B^2+b3B''2+b4B^2; 

b5C=blC'^2+b2C''2+b3C'^2+b4C'^2; 

b5D=blD^2+b2D^2+b3D^2+b4D^2; 

b5E=blE^2+b2E'^2+b3E^2+b4E''2; 

b5F=blF'^2+b2F'^2+b3F''2+b4F'^2; 

b5G=blG'^2+b2G'^2+b3G''2+b4G'^2; 

b5H=blH'^2+b2H'^2+b3H''2+b4H'^2; 

b5I=blH''2+b2H'^2+b3H^2+b4I^2; 

b5J=blJ'^2+b2J^2+b3J''2+b4J'^2; 

b5K=blKA2+b2K''2+b3K'^2+b4K'^2; 

b5L=blL'^2+b2L'^2+b3L'^2+b4L'^2; 

b5M=blM''2+b2M^2+b3M'^2+b4M^2; 

b5N=blN'^2+b2N'^2+b3N''2+b4N'^2; 

b50bl0^2+b2C>'2+b30^2+b40'2; 

b5P=blP^2+b2P''2+b3P''2+b4P^2; 


b00=.2182;  %  approximate  zero 
bl0=.9759;  %  approximate  1 

%  Determine  k  constant  for  the  variance  of  X 

nu=.l; 

delk=.3; 

ka=(2*pi*nu)/(4*pi''2*(nu''2+delk^2)); 

kb=l/(16*piM*(nu''2+<ielk''2y'2); 

kc=(2*pi*nu*exp(-2*pi*nu))*(2*pi*deIk*sin(2*pi*delk)-2*pi*nu*cos(2*pi*delk)); 

kd=(2*pi*delk*exp(-2*pi*nu))*(2*pi*(delk*cos(2*pi*delk)+nu*sin(2*pi*delk))); 

ke=4*pi''2*(nu''2-delk''2); 

k=ka-kb*(kc-t-kd+ke); 

n*=l; 

j=i; 

m=l; 

%  -60dB<=Peak  Power<=-44dB 
%  Peak  power  is  A^2/2  =  Ps 
for  ss=-60:.5:-44; 
ps(m^  Ky'Css/lO); 
m=m+l; 
end 


77 


%  The  goal  here  is  to  solve  the  equation  for  the  optimal  threshold 

%  vs  power 

alpha>Bzeros(length(ps)Jength(ck)); 

PBEBzeros(len^(ps)Jen^(ck)); 

for  j«l :length(ck);  %  coupling  values  loop 
n^l; 

for  n«l  :length(ps);  %  peak  ppower  values 
Ps«ps(n); 

%  Constants 

qcns^q/(Ps*T*R); 

vnoise=N0/(R''2*Ps^2*T); 

v(iark*(q*Idk)/(R'^2*Ps'^2*T); 

%  Determine  the  sigma  for  each  of  the  0  bit  patterns 

s0A=S(pt(4*ck(j)*b00^2*k*b5A+qcnst*(b0(y'24ck(j)*b5A)+vnoisc+vdaik); 

s0B=sqrt(4*ck0)*b00^2*k*b5B4qcnst*(b00''2+ck0)*b5B)+vnoise+vdaTk); 

s0C=sqn(4*ck(j)*b00^2*k*b5C+qcnst*(b00^2+ck0)*b5C)+vnoise+vdark); 

sOD=sqrt(4*ck0*bOO^2*k*b5D+qcnst*(bOO''2+ck(])*b5D)+vnoise+vdark); 

s0E*sqrt(4*ck(j)*b0(y'2*k*b5E4qcnst*(b00^2+ck(j)*b5E)+vnoise+vdark); 

sOF=sqrt(4*ck^)*b(XK'2*k*b5F+qcnst*(bOO'^2+ckO)*b5F)+vnoise+vdarkV, 

sOG=sqrt(4*ck0*b(Xy^2*k*b5G4qcnst*(b(Xy^2+ck(j)*b5G)+vnoise+vdark); 

s0H=sqrt(4*ckQ)*b0(y^2*k*b5H+qcnst*(b00'^2+ck0)*b5H)+vnoise+vdark); 

s0I=sqrt(4*ck(j)*b00^2''’k*b5I+qcnst*(b0(y^2+ck(j)*b51)+vnoise+vdaA); 

s0J=sqrt(4*ck0‘)*b00^2*k*b5J+qcnst*(b0(y^2+ck0)*b5J)+vnoisc+vdark); 

s0K=sqrt(4*ck(j)*b00^2*k*b5K+qcnst*(b00^2+ck(j)*b5K)+vnoise+vdark); 

sOL=sqrt(4*ck(j)*bO(y'2*k*b5L+qcnst*(b(Xy'2+ck(j)*b5L)+vnoise+vdark); 

sOM=sqrt(4*ck(j)*bO(y^2*k*b5M+qcns-'*‘(bOO'^2+ck(j)*b5M)+vnoise+vdaA); 

sON=sqrt(4*ck(j)*bO(y'2*k*b5N+qcnst*(b()(y^2+ck(j)*b5N)+vnoise+vdark); 

s0O=s<pt(4*ck(j)*b00^2*k*b5C)+qcnst*(b0(y'2+ck(j)*b5O)+vnoise+vdark); 

s0P=sqrt(4*ck(j)*b(XK'2*k*b5P+qcnst*(b00^2+ck(j)*b5P)+vnoise+vdark); 

%  Determine  the  variance  for  each  of  the  0  bit  patterns 

sOA2=sOA''2; 

sOB2=sOB''2; 

sOD2=sOD^2; 

sOE2=sOE''2; 

sOF2*sOF''2; 

sOG2ssOG^2; 

sOH2=sOH'^2; 

sOI2=sOI-'2; 

sOJ2=sOJ^2; 

sOK2-sOK''2; 

sOL2»sOL''2; 

sOM2=sOM''2; 

sON2=sON'^2; 

s002=s00^2; 

sOP2=sOP''2; 
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%  Determine  the  sigma  for  each  of  the  1  bit  patterns 

slA«sqrt(4*ck(j)*bl(y^2*k*b5A+qcnst*(bl0^24ck(j)*b5A)+vnoise+vdark); 

slB=sqr^4*ckO)*bl(y^2*k*b5B+qcnst*(bl(y'2-K:k(j)*b5B)+vnoisc+vdaik); 

slC»sqrt(4*ck(])*bl(y'2*k*b5C4qcnst*(bl(y'2-fck(j)*b5C)+vnoisc+vdaik); 

slD»s<pt(4*ckQ)*bl0^2*k*b5D+qcnst*(bl0^2+ck(j)*b5D)+vnoisc+vdaik); 

slE*sqrt(4*ckO)*bl(y^2*k*b5E+qcnst*(bl0^2+ck(j)*b5E>+vnoise+vdark); 

slF*sqrt(4*ck^)*bl(y'2*k*b5F+qcnst*(bl0^2+ckO)*b5F)+vnoisc+vdark); 

slG=sqn(4*ck(j)*bl(y'2*k*b5G-Kjcnst*(bl(y'2+ck(j)*b5G)+vnoise+vdark); 

slH»sqn(4*ck(j)*bl0^2*k*b5H4qcnst*(bl0^2+ck(j)*b5H)+vnoisc+vdark); 

slI=s(pt(4*ck(j)*bl0^2*k*b5I+qcnst*(bl0^2+ck(j)*b5I)+vnoise+vdaik); 

slJ*=sqrt(4*ck(j)*bl(y'2*k*b5J+qcnst*(bl(y^2+ck(j)*b5J)+vnoisc+v(iark); 

slK=sqrt(4*ck(j)*bl0^2*k*b5K+qcnst*(bl0^24ck(j)*b5K)+vnoise+vdark); 

slLsssqrt(4*ck(j)*bl(y^2*k*b5L+qcnst*(bl0^24ck(j)*b5L)+vnoisc+vdark); 

slM=sqrt(4*ck(j)*bl0^2*k*b5M+<}cnst*(bl0^24ck(j)*b5M)+vnoise+vdark) 

slN=sqrt(4*ck(j)*bl0^2*k*b5N+qcnst*(bl(y'2+ck(j)*b5N)+vnoisc+vdark); 

slO=sqrt(4*ckO)*bl0^2*k*b5CHqcnst*(blO''2+ck(j)*b50)+vnoise+v(iaik); 

slP=s<frt(4*ck(j)*bl0^2*k*b5P+qcnst*(bl0^2+ck(j)*b5P)+vnoise+vdark); 

%  Determine  the  variance  fw  each  of  the  1  bit  patterns 

slA2=slA'^2; 

slB2=slB''2; 

slC2=slC''2; 

slD2=slDA2; 

s1E2=s1EA2; 

slF2=slF''2; 

slG2*slG''2; 

slH2=slH''2; 

sll2*sll'^2; 

slJ2=slJ''2; 

slK2=slK^2; 

slL2=slL'^2; 

slM2=slM''2; 

slN2«slN''2; 

sl02=slO2; 

slP2=slP'^2; 


%  Determine  the  mean  for  each  of  the  0  bit  patterns 

mOA=b00^24ck(j)*b5A; 

m0B=b00^2+ck(j)*b5B; 

mOObO(y'2+ck(j)*b5C; 

m0D=b00^24ck(j)*b5D; 

mOE=bOO^24ck0)*b5E; 

mOF=bOO^24ck0‘)*b5F; 

m0G=b00^2+ck(j)*b5G; 

mOI=bOO^2+ck0)*b5I; 

m0J=b00^2+ck(j)*b5J; 

m0K=b00^2+ck(j)*b5K; 

mOL=b00^2+ck(j)*b5L; 

m0M»b00^2+ck(j)*b5Mi 

m0N*b00^24ck(j)*b5N; 
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niOO-bO(y'24ck(j)*b5a 

m0P«b00^2-fck(j)*b5P; 


%  Detmnine  the  mean  for  each  of  the  1  bit  patterns 

mlA=bl0^2+ck(j)*b5A; 

mlB»bl0^2+ckO)*b5B; 

mlC=blOA2+ck(j)*b5C; 

mlD=bl0^24ck(j)*b5D; 

mlE=blO^2+ck0')*b5E; 

mlF«bl0^24ck(j)*b5F: 

mlG=bl0^24ck(j)*b5G; 

mlH=bl0^24ck(j)*b5H; 

mlI«bl0^2-K:k(j)*b5I; 

mlJ=bl0^2+ck(j)*b5J; 

mlK=bl0^24ck(j)*b5K; 

mlL=bl0^2+ck(j)*b5L; 

mlM=bl0^2+ck(j)*b5M; 

mlN=bl0^2+ck(j)*b5N; 

ml0bl0^2+ck(j)*b50; 

mlP=bl0^24ck(j)*b5P; 

%Solve  the  equation  F(a)=0  to  optimize  a,  the  threshold.  F(a)aO  comes 
%  from  (dP/da)=0 

a=0:.0009:l; 

FA=(l/s0A)*exp(-((a-m0A).''2)./(2*s0A2))-(l/slA)*exp(.((mlA-a).^2)/(2*slA2)); 

FB=(l/s0B)*exp(-((a-m0B).'^2)./(2*s0B2))-(l/slB)*exp(-((mlB-a).^2)/(2*slB2)); 

FC=(l/s0C)*exp(-((a-m0C).''2)./(2*s0C2))-(l/slC)*exp(.((mlC-a).'^2)/(2*slC2)); 

FD=(l/s0D)*exp(-((a-m0D).'^2)./(2*s0D2))-(l/slD)*exp(-((mlD-a).^2)/(2*slD2)); 

FE=(l/s0E)*exp(-((a-m0E)/2)./(2*s0E2)Hl/slE)*exp(-((mlE-a).^2y(2*slE2)); 

FF=(l/s0F)*exp(-((a-m0F).^2)./(2*s0F2))-(l/slF)*exp(-((mlF-a).''2)/(2*slF2)); 

FG=(l/s0G)*exp(-((a-m0G).''2)./(2*s0G2))-(l/slG)*exp(-((mlG-a).^2)/(2*slG2)); 

FH=(l/s0H)*exp(-((a-m0H).''2)./(2*s0H2))-(l/slH)*exp(-((mlH.a).'^2)/(2*slH2)); 

n=(l/s0I)*exp(-((a-m0I).''2)./(2*s0I2))-(l/slI)*exp(-((mlI-a).^2)/(2*slI2)); 

FJ=(l/s0J)*exp(-((a-m0J).''2)./(2*s0J2))-(l/slJ)*exp(-((mlJ-a).''2)/(2*slJ2)); 

FK=(l/s0K)*exp(-((a-m0K).^2)./(2*s0K2))-(l/slK)*exp(-((mlK.a).^2)/(2*slK2)); 

FL=(l/s0L)*exp(-((a-m0L).''2)./(2*s0L2))-(l/slL)*exp(-((mlL-a).'^2)/(2*slL2)); 

FM=(l/sOM)*exp(-((a-mOM).''2)y(2*sOM2))-(l/slM)*exp(-((mlM- 

a).''2)/(2*slM2)); 

FN=(l/s0N)*exp(-((a.m0N).^2)./(2*s0N2))-(l/slN)*exp(-((mlN-a).'^2)/(2*slN2)); 

F0=(l/s00)*exp(-((a-m00).''2)./(2*s002))-(l/sl0)*exp(-((miaa).''2)/(2*sl02)); 

FP*(l/s0P)*exp(-((a-m0P).'^2)7(2*s0P2)Hl/slP)*exp(-((mlP-a).^2)/(2*slP2)); 

F=FA+FB+FC+FD+FE+FF+FG+FH+n+FJ+FK+FL+FM+FN+FO+FP; 

[Y,I]*min(abs(F)); 


alpha(nj)=a(I); 

aa=alpha(nj); 

%  Determine  &e  probability  of  bit  error  for  each  optimal  threshold 
%  over  all  bit  panems  A,BfCJ) 
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pOA(n)>«rf((alpha(iiJ)-niOA)/sOA4nO; 

pOB(n)Berf((alpha(nj)-inOBVsOB3nf); 

pOC(n>Berf((alpha(nj)-mOC)/sOC4nf); 

pOD(n>>erf((alpha(nJ)-inOD)/sOD4nO: 

pQE(n)aerf((alpha(nj)-inOE)/sOE4nf): 

pOF(n)=crf((alpha(nj)-inOFysOF4nf); 

pOG(n)seif((alpha(nJ)-inOG)/sOG4nf): 

pOH(n)a^((alpha(nJ)-inOH)/sOH4nf): 

pOI(n)«erf((alpha(nj>mOI)/sOI4nf): 

pOJ(n)«af((alpha(nJ)-mO^/sOJ4nO; 

pOK(n)aeif((ajpha(nJ)-nk)K)/sOKinf); 

pOL(n)aaf((alpha(nj)-mOL)/sOL4nf): 

pOM(n)setf((a!idia(nJ)-mOMVsOM4^: 

pON(n)serf((alpha(nJ)-mON)/sON,inf); 

p00(n>=erf((alpha(nj)-in00)/s004nf); 

pOP(n)=crf((alpha(nJ)-inOPi/sOP,inf); 

p  1  A(n)=crf((ml  A-alpha(n  j))/s  1  A, inf); 

plB(n)serf((mlB-alpha(nJ))/slB,inf); 

plC(n>=crf((mlC-alpha(nj)yslC4nf); 

plD(n)=erf((inlD-alpha(nJ)VslD,inf); 

plE(n)=crf((inlE-alphaCnj))/slE4nO; 

p  lF(n)=erf((m  1  F-alpha(n  j))/s  1  F4nf); 

plG(n)=«if((mlG  alpha(n  j))/s  lG4rtO; 

plH(n)=cif((mlH-alpha(nj)yslH4nf); 

pll(n)serf((mll-alpha(nj))/sll4nf); 

plJ(n)=«rf((inlJ-alpha(nj))/slJ4nO; 

plK(n)=cif((inlK-^pha(no))/slK4nO; 

plL(n)=eif((mlL-alpha(nj))/slL4nf); 

plM(n)=serf((nilM-ipha(nj))/slM4n]f); 

plN(n)=^((inlN-alpha(no))^lN4nf); 

plO(n)=erf((mlO-alpha(nj))/slO,inO; 

plP(n)=crf((mlP-alpha(nJ))^lP4nf); 


pA(n)=pOA(n)+plA(n); 

pB(n)^)OB(n)+p  1  B(n); 

pC(n)^)OC(n)+plC(n); 

pD(n)=4>0D(n)+plD(n); 

pE(n)=pOE(n)+plE(n); 

pF(n)=pOF(n)+p  lF(n); 

pG(n)=pOG(n)+p  1  G(n); 

pH(n)=pOH(n)+plH(n); 

pl(n)=j01(n)+pll(n); 

pJ(n)=pOJ(n)+plJ(n); 

pK(n)=pOK(n)+plK(n); 

pL(n)*^L(n)+plL(n); 

pM(n)^M(n)+plM(n); 

pO(n)=pOO(n)+plO(n); 

pP(n)=pOP(n)+plP(n); 

PBE(nj)=(pA(n)+pB(n)+pC(n)+pD(n)+pE(n)+pF(n)+pG(n)+pH(n)+pI(n)+pJ(n) 

+pK(nHpL(n)+pM(n)+pN(n)+pO(n)+pP(n))/64; 
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end  %Psloop 
end  %Qckxq) 


al>>«lpha(:,l); 

a2Balpha(:t2y, 

a3«alpha(:3); 

a4«alpha(:,4); 

a5>alpha(:,5); 

a6^pha(:,6); 

a7asalpha(:,7); 

a8=a]pha(:,8); 

a9Balpha(:.9); 

alO<Balpha0.1O); 

all>idpha(:,ll); 

pl«(PBE(;.l)); 

p2=(PBE(:.2)); 

p3=(PBE(:.3)); 

p4=(PBE(:,4)); 

p5=(PBE(:.5)); 

p6=(PBE(:,6)); 

p7=(PBE(:.7)); 

p8=(PBE(:.8)); 

p9=(PBE(:.9)); 

pl0=PBE(:,10); 

pll=PBE(:,ll); 


PB=10*logl0(ps); 

plot(PB,al,'-’,PB.a2;-’.PB.a3/-’,PB,a4  v,PB,a5.’-’.PB,a6,' 
grid... 

%title('case  1.3,nu=,l,delk=.3:  Four  channels'),.. 
xlabeK'Peak  Powcr(dBW)’),ylabel('Optimal  Threshold') 


’.PB,a7,'- 


s=['Crosstalk=’]; 

texi(.15,.6,s,'sc') 

cl=0; 

s=['  ',num2str(cl)]; 
text(.2,.55,s,'sc') 
[xs,ys]=dc2sc(PB  (19),al(19)); 
’'olyline([.3,xs],[.56,ys],'-r','sc') 

c2=-35; 

s=[nuni2str(c2),'  dB']; 
text(.2,.52,s,*sc') 
[xs,ys]=dc2sc(PB(19),a2(19)); 
polyline([.3,xs],[.53,ys],'-r’,'sc') 


c3=-32; 

ss[nui:]^str(c3),'  dB']; 
text(.2,.49,s,'sc’) 
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[xs.ysJ-dc2sc(PB(19).a3(19)); 

polylinc([.3,xs],[.50,ys],'-r',’sc') 

c4=-29; 

s«[nuiD‘>str(c4),'  dB’]; 
text(.2..4^,s,'sc’) 
[xs.ys]«<ic2sc(PB(19).a4(19)); 
polyline([.3,xs],[.47,ys],'-r*.'sc') 

c5^26; 

ss[nuin2str(cS),'  dB']; 
tcxt(.2,.43,s,'sc') 
[xs.ys]=dc2sc(PB(19).a5(19)); 
polyline([.3,xs],[.44,ys],'-r’,'sc*) 

c6=-23; 

ss[nuin2str(c6),'  dB']; 
tcxt(.2,.40,s,'sc') 
[xs,ys]=dc2sc(PB(19).a6(19)); 
polyline([.3.xs],[.41,ys].'-r','sc') 


c7=-20; 

s=[num2str(c7),'  dB']; 
text(.2,.37,s,'sc') 
[xs,ys]=dc2sc(PB(19).a7(19)); 
polyline([.3,xs],[.38,ys],'*r','sc') 


meta  t4chl3xa 
pause 


axis([>60  -44  -21  -5]) 

plot(PB,pl,’-'.PB,p2.'-',PB.p3.'-',PB,p4.'-'.PB,p5,'-',PB.p6.' 

%PB.p8,'-',PB,p9.'-',PB,plO,'-'.PB.pll.'-') 

scmilogy.grid,.. 

%title('case  1.3,nu=.l,dclk=.3:  Four  channels'),.. 
xlabelCPeak  Power(dBW)'),ylabel('Probability  of  Bit  Error') 


',PB,p7.'-')... 


s=['Crosstalk=']; 

text(.15,.6,s,'sc') 


cl*-20; 

s=[nuin2str(cl),'  dB']; 
text(.2,.55,s,'sc') 
[xs,ys]=dc2sc(PB(15).p7(15)); 
polyline([.3,xs],[.56,ys],’-r',’sc') 

c2=-23; 

s=[num2str(c2),'  dB']; 
text(.2,.52,s,'sc’) 
[xs,ys]*dc2sc(PB(15),p6(15)); 
polyline([.3,xs],[.53,ys],'-r','sc') 

c3*-26; 
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s«[num2str(c3),'  dB']; 
tcxt(.2,.49,s,'sc’) 
[xs,ys]-dc2sc(PB(15).p5(15)); 
polyline([.3,xsl,(.50,ysl,’-r’,*sc') 

c4a^29; 

s«[nua^str(c4),'  dB']; 
tcxt(.2,.46,s,'sc’) 
[xs.ys]-dc2sc(PB(15).p4(15)); 
polylinc([.3,xs],I.47,ys],'-r‘,’sc') 

c5*-32' 

s>:[num2str(c5),'  dB']; 
text(.2,.43.s,'sc') 
[xs,ys]»dc2sc(PB(15).p3(15)); 
polylinc([.3,xs],[.44,ys],'-r','sc') 

c6=-35; 

ss[nuin2str(c6),'  dB']; 
tcxt(.2,.40,s,'sc') 
[xs,ys]»dc2sc(PB(15).p2(15)); 
polylinc([.3,xs],[.41,ys],'-r','sc') 

c7=0; 

s=['  'jjum2str(c7)]; 
text(.2,.37,s,'sc') 
[xs.ys]=dc2sc(PB(  1 5).pl  ( 1 5)); 
polylinc([.3,xs],[.38,ys],'-r','sc') 
axis; 

meta  t4chl3xb 
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